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ABSTRACT 

A  detailed  study  is  presented  of  the  acoustic  radiation 
pressure  exerted  by  plane  compress ional  waves  in  non- viscous 
liquids  and  gases  upon  a  plane  obstacle.  The  present  report  is 
largely  a  further  development  and  extension  of  a  very  comprehensive 
and  penetrating  treatment  of  acoustic  radiation  pressure  by  L« 
Brillouin.  The  theory  of  plane  waves  in  liquids  is  extended  to 
the  case  of  finite  amplitudes.  From  this  more  general  point  of 
view  the  effects  at  small  amplitudes  are  derived  and  discussed  in 
detail.  Formulas  are  given  for  normal  and  oblique  incidence  of 
the  acoustic  bean  for  small  amplitudes s  valid  for  any  reflection 
coefficient  of  the  receiving  plane.  The  radiation  pressure  at 
finite  amplitudes  upon  a  perfect  absorber  is  calculated. 

Special  consideration  is  devoted  to  the  actual  physical 
processes  involved;  the  meaning  of  the  so-called  Rayleigh  pressure 
and  Langevin  pressure  is  discussed.  For  gases  the  radiation 
pressure  as  well  as  the  Rayleigh  pressure  in  progressive  and  standing 
waves  are  computed. 

The  report  concludes  with  an  assembly  of  the  chief 

equations. 


ICC  |pqft->|iq  «n  4k!  e=^ 


SYMBOLS 
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a  =  space  coordinate 

C  =  constant  of  integration 


P 


E 


in 

ot 


phase  velocity  of  light  or  sound 
specific  heat  at  constant  pressure 

tine-average  value  of  the  density  of  total  energy  for  a 
harmonic  v/ave  traveling  in  one  direction  (E  =  poco^Q^/2) 

tine-average  of  the  density  of  acoustic  kinetic  energy 

tine-average  of  the  density  of  acoustic  potential  energy 
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k 

l  - 

El  = 


Prads  P 


x 


p 

pt 

S 

T 


T 

lk 


u  su  9U 
Xs  ys  2 


electric  field  strength 
frequency 

magnetic  field  strength 

width  of  surface  receiving  acoustic  radiation  (Fig.  12) 

intensity  (power  per  unit  area)  of  acoustic  radiation  per 
square  centimeter 

Bessel  function  of  order  1  and  argument  x 
=  oo/c 

length  of  surface  receiving  acoustic  radiation  (Fig,  12) 
mass 

=  radiation  pressure  or  its  component  in  the  positive 
x-direction 

hydrostatic  pressure 

excess  pressure  due  to  compress ional  wave  motion 
pQ  +  p  =  total  pressure 

shearing  stress  due  to  radiation  pressure  at  oblique  incidence 

stress  tensor,  also  absolute  temperature 

components  of  stress  tensor 

l/f  =  period 

tine  coordinate 

=  particle  velocities  in  the  x-?y-sz -directions 


u 


x»y,z 

xSy* 


•  *  • 
111 

particle  velocity  in  a  plane  conpressional  nave  along 

the  x-direction  (u  =  u  ) 

x 

space  coordinates 

space  coordinates  of  a  system  rotated  through  an  angle  gfi. 

(Fig.  5) 


V  =  unit  volume;  ^  V  =  snail  change  in  V 

VQ  =  undisturbed  unit  volume 
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T° 
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/,v- 


(0 
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cr 


compressibility  of  liquid  (cm^  per  unit  force) 

amplitude  reflection  coefficient  ( ^  =  1 ) 

cp/cv  =  ratio  of  specific  heats  at  constant  pressure  (c  ) 
and  constant  volume  cy) 

(c  +  c  +  e  ) 
xx  yy  zz 

amplitude  transmission  coefficient  ( <f  -  1 ) 
dielectric  permittivity 

components  of  strain  tensor  (see  Eq.  (2.7)) 
viscosity  coefficient 


=  phase  angle  of  reflected  wave 

=  angle  between  the  normal  to  a  receiving  plane  and  the 
direction  x  of  the  \;ave  propagation 

=  wavelength 

=  magnetic  permeability 

=  density;  jfQ  =  undisturbed  density 

=  entropy 


Tty  =  cot  +  lex 

rj9  C  =  particle  displacements  in  the  x-  ,y-,,z- directions 
cn  =  2irf  =  angular  frequency 


Ilote 

One  bar  over  a  symbol  denotes  the  time-average  value  of  the 
quantity  concerned,  as  p,  two  bars  the  average  value  in  time  and 
and  spaces  as  p. 


A  star  *  indicates  that  the  quantity  concerned  refers  to  a 
moving  particle  or  volume  element  (Lagrangian  coordinates),  as  p*» 
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PART  I 

THE  CCDCEFT  OP  RADIATIOR  PRESSURE 

1  *  Radiation  Pressure  in  Electrodynamics 
The  concept  of  radiation  pressure  originated  in  electro¬ 
dynamics.  If  we  consider  a  plane  surface  emitting  a  plane 
electromagnetic  v/ave  (Fig.  l) 
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Fig.  1 

Plane  electromagnetic  wave  emitted 
by  a  plane  surface 

in  the  positive  redirection,  a  reacting  force  is  exerted  upon  the 
emitting  plane  due  to  the  transport  of  momentum  by  the  electro¬ 
magnetic  field  in  the  direction  of  nave  propagation.  Imagine  a 
cylinder  with  a  cross-section  of  1  cm  perpendicular  to  the  x-axis, 
A  nave  front  leaving  the  emitting  surface  at  any  time  t  reaches 
a  cross-section  of  this  cylinder  at  the  distance  c  one  second  later, 
where  c  is  the  velocity  of  the  wave.  Assuming,  for  simplicity, 
an  electromagnetic  wave  of  a  rather  high  frequency,  so  great  a 
number  of  wavelengths  may  fill  the  distance  c  that  we  can  speak  of 
a  mean  total  electromagnetic  energy  averaged  over  the  length  c  of 
our  cylinder.  Dividing  this  mean  total  energy  by  c  we  obtain  the 
mean  energy-density  that  is,  the  mean  energy  per  cm^  of  the 
electromagnetic  wave.  If  the  emission  is  steady,  E  is  independent 


of  time. 
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The  whole  energy  filling  the  cylinder  of  length  c,  and 
therefore  leaving  the  emitting  surface  each  second.,  is  E*c,  Accord¬ 
ing  to  the  equivalence  of  mass  and  energy  this  corresponds  to  a  mass 
transport  in  the  positive  x-direction.  Let  the  whole  mass  filling 
the  cylinder  of  length  c  be  a  =  po,  p  denoting  the  electromagnetic 
"mass-density1*  per  cmr .  Then  we  have 


and 


(1.1) 


The  equivalent  mechanical  momentum  emitted  each  second  from  each 
cm^  of  the  radiating  surface  is  m*c  =  (pc)*c  =pc^  =  E.  Accord¬ 
ing  to  Newton’s  Law  a  reaction-force  is  exerted  on  each  square 
centimeter  of  the  radiating  surface,  which  equals  the  momentum 
me  =  E  contributed  per  second  to  the  emitted  wave.  This  reacting 
force  Px  per  cn^  in  the  negative  x-direction  is  numerically  equal 
to  the  so-called  electromagnetic  radiation  pressure  Prad: 


Px  me  =  Ppad  E 


(1.2) 


Prad  equals  the  mean  electromagnetic  density  averaged  over 
a  wavelength  in  space: 

Pr„d  =  E  =  ~  (e  -  +  M  -  )  dx  (1.3) 

A  da  2  2 

where  X=  c/f  =  wavelength;  c  is  the  permittivity;  jj.  is  the 
permeability;  IS,  if  are  the  electric  and  magnetic  field  vectors 


of  the  plane  wave. 
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If  the  emitting  surface  were  free  to  move,  this  radiation- 
force  would  cause  a  notion  of  the  surface  in  the  direction  opposite 
to  that  of  the  emitted  radiation.,  in  such  a  way  that  the  mass- 
center  of  the  whole  system  (emitting  mass  +  mass -equivalent  of 
radiated  energy)  remained  fixed  in  space. 

How  we  consider  the  case  in  which  the  emitted  radiation 
strikes  a  plane  totally  absorbing  surface  perpendicular  to  the  x- 
axis.  Such  a  surface  is  usually  called  a  surface  of  a  black  body. 
The  electromagnetic  momentum  which  one  cn^  of  the  black  surface 
absorbs  each  second  amounts  likewise  to  me  ~ j&c  ~  E.  This  gain 
in  momentum  per  second  corresponds  to  a  radiation  force  on  each 
cm^  of  the  receiving  area  in  the  positive  x-direction.  So  xje  say 
that  the  electromagnetic  radiation  exerts  a  radiation  pressure  upon 
each  square  centimeter  of  the  absorbing  surface,  which  equals 
the  mean  energy  density  E  of  the  electromagnetic  wave. 

If  the  surface  perpendicular  to  the  direction  of  motion 
of  the  electromagnetic  wave  is  not  black,  but  totally  reflecting, 
the  momentum  of  the  wave  changes  its  sign  at  the  reflector  from 
+  me  to  -  me.  So  the  total  change  in  momentum  per  second  and  per 
cm^  of  the  wave  amounts  to  2nc  =  2  pc^  =  2E.  Therefore  the  radia¬ 
tion  pressure  Pracj  on  a  perfectly  reflecting  surface  is  equal  to 
2E.  If  the  receiving  surface  partially  absorbs  and  partially  re¬ 
flects,  the  radiation  pressure  has  a  value  between  £  and  2E, 
depending  on  the  coefficient  of  the  reflection  of  the  surface, 
which  expresses  the  ratio  of  the  reflected  to  the  incident  x/ave 
energy. 

2,  The  Concent  of  Pressure  in  Fluids 

In  view  of  our  later  treatment  of  acoustic  radiation 


pressure*  we  first  consider  the  meaning  of  the  word  npressureu. 

In  the  physics  of  deformable  bodies  (solids*  liquids*  gases) 
the  word  pressure  is  exemplified  by  the  hydrodynamic  pressure  in 
a  fluid.  The  pressure  is  a  scalar  p  (x,y*z )  which  nay  change 
from  point  to  point  in  space.  If  we  consider  a  volume-element 
of  arbitrary  shape  inside  a  fluid*  the  hydrodynamic  pressure  p 
exerts  the  sane  force  in  all  directions  if  tlie  volume  is  so  small 
that  the  change  in  pressure  with  respect  to  x*y*z  can  be  neg¬ 
lected  (Fig.  2)-, 


Fig.  2 

Hydrostatic  pressure  p  in  fluids 

Or*  if  we  imagine  a  small  surface  element  in  the  fluid*  the  pressure 
p  on  it  is  the  same  for  each  orientation  in  space  and  always  normal 
to  the  surface. 

Instead  of  using  the  word  pressure  we  can  also  say  that 
there  is  a  stress  T  acting  on  any  small  element  of  area  perpendic¬ 
ular  to  the  element  and  of  the  same  value  for  every  orientation  in 
space  of  the  element.  fie  write 

T  =  -  p  (1.4) 

since  usually  the  pressure  is  called  positive  if  a  volume  element 
is  compressed*  while  a  stress  is  called  positive  when  tensile 
(Fig.  3). 
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Pig.  3 

Directions  of  positive  pressure  p 
and  positive  uniform  stress  T 


At  this  point  it  is  desirable  to  consider  the  general  ex¬ 
pression  for  a  stress-tensor.  This  expression  will  then  be 
specialized  for  the  case  of  a  fluid.  Any  stress,  whether  due  to 
mechanical  forces  or  to  an  electromagnetic  field,  is  a  tensor  hav¬ 
ing  in  general  nine  components.  It  is  represented  by 


xik  “1  1 
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XX 

xy 

xz 
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Tyy 
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zx 

zy 

zz 

(1.5) 


The  T..  (T  T  .  T  )  represent  the  normal  forces  perpendicular 
ii  xx  yy  zz  ■ 

to  surface  elements  in  the  yz-,  xz-,  and  xy-p lanes ,  whereas  the 
shearing  forces  T.^  (T^,  •••)  act  parallel  to  these  surface 

elements.  In  fluids  under  mechanical  action,  disregarding  viscous 
forces,  only  the  T..  exist  and  the  stress  tensor  reduces  to 
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XX 
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0 

r?  — 
"liq 
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yy 

0 
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(1.6) 


Furthermore,  as  we  have  stated,  the  force  on  a  surface  element  in  a 
fluid  is  independent  of  its  orientation  in  space,  so  that 
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or 
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(1.7) 


According  to  the  rules  of  tensor  calculus  the  quantity 


P  +  T  +  T 
:oc  yy  zz 


(1.8) 


is  invariant 5  that  is,  it  has  the  same  value  for  any  set  of  rec¬ 
tangular  coordinates  at  the  sane  point.  In  our  case  T„r  =  T™  ~ 

Jwv  j  J 

TZz  =  -  p,  and  the  quantity  I  =  -  3p,  or 


-  P  = 


2 


(1.9) 


The  hydrostatic  mean  pressure  p,  which  is  of  course  independent  of 
the  system  of  coordinates  chosen,  can  therefore  in  general  be 
defined  by 

(TSC  +  Tjrjr  +  '  (1-10) 

3  •  The  Tensorial  Character  of  Electrodynaraic  Radiation 
Pressure 

In  considering  electrodynaraic  radiation  pressure,  we  found 
that  the  force  P  due  to  a  wave  propagated  in  the  x-direction  is 

.X. 

also  directed  in  the  x-direction*  Ho  force  is  acting  in  any 
direction  perpendicular  to  the  x-axis.  So  the  radiation  pressure 


is  not  a  "pressure"  in  the  sense  in  which  we  use  this  notion  in 
hydrodynamics.  Indeed,  if  we  change  the  direction  of  our  re¬ 
ceiving  surface  in  such  a  way  that  its  outer  normal  makes  an 


angle  S  with  the  positive  x-axis  (Fig.  4)s 


\  X-*1 

_ 1  \,  -^3- 

' .  'f\  x 
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/ 

Fig.  4 

Oblique  incidence 


the  radiation  force  perpendicular  to  a  perfectly  reflecting  surface 
is  known  to  be 

Px»  =  2E  cos2^  (1.11) 

It  is  by  no  neans  independent  of  the  orientation  of  the  surface  with 
respect  to  the  incident  wave.  This  neans  that,  strictly,  one 
should  speak  of  the  radiation  tensor,  which  at  normal  incidence  has 
only  one  component,  =  -  Px.  As  only  energy-densities  are  in¬ 
volved  in  p  =  E,  the  expression  for  T  is  independent  of  the  polar¬ 
ization  of  the  plane  electromagnetic  wave: 


(1.12) 


According  to  the  rules  for  the  transformation  of  a  tensor  T  to  a 


new  set  of  axes,  whose  x,  y-coordinates  are  rotated  through  an 
angle  d?  about  the  z-axis  (?ig.  5)j  the  tensor 
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But  according  to  Eg,  (1.1 3),  the  radiation  tensor  (1.1 2)  transforms 
thus : 


p  2  A 

cos  2K 

Y  sin 

0 

P 

,  - 

T’rad  = 

Y  sin  2  2Jp 

| 

-Px  sin^z$ 

0 

1 

o 

1 

0 

Cj 

Ue  have  discussed  the  difference  between  the  concept  of 
a  hydrodynamic  pressure,  which  is  the  same  in  all  directions,  and 
the  physical  properties  of  the  electrodynamic  radiation  pressure, 
in  some  detail,  because  it  will  turn  out  that  also  in  the  case  of 
acoustic  radiation  pressure  the  tensorial  character  of  this  quan¬ 
tity  must  be  taken  into  account. 

The  electromagnetic  radiation  pressure  is  proportional 
to  the  mean  total  energy  density  and  therefore  to  the  square  of 
the  amplitudes  of  the  electric  and  magnetic  field-strengths 
(Eq.  1.3).  This  holds  for  all  field-strengths,  whether  small  or 
large.  As  to  the  acoustic  radiation  pressure,  we  shall  see  that 
a  similar  law-  is  valid  only  for  sufficiently  small  amplitudes  of 
vibration. 

4o  Radiation  Pressure  in  Acoustics 

Rayleigh  was  the  first  to  apply  the  concept  of  radiation 
pressure  to  mechanical  waves  in  gases. ^  He  established  relations 
for  the  average  value  in  time  of  the  pressure  produced  by ‘an  infinite 
plane  wave  in  fluids  and  showed  that  this  mean  pressure  is  pro¬ 
portional  to  the  mean  energy  density  of  the  wave  motion.  The 
factor  of  proport  tonality  was  found  by  him  not  to  be  one  in  general, 
as  in  the  case  of  electromagnetic  waves,  but  dependent  on  the 


14) 


TO 


special  lav;  connecting  the  pressure  p  and  the  density  p  in  the  fluid 
under  consideration.  The  physical  quantity  called  by  Rayleigh 
the  ’’pressure  of  vibrations",  is,  however,  not  identical  with  what 
is  measured  usually  as  "radiation  pressure".  The  "radiation 
pressure"  of  a  compressional  wave  striking  a  perfect  absorber  equals 
exactly  twice  the  mean  kinetic  energy  density  of  the  wave  motion, 
as  we  shall  see.  For  small  amplitudes  this  expression  becomes 
identical  with  the  mean  total  energy  density  and  we  obtain  in  this 
case  the  sane  relation  as  in  the  electromagnetic  case. 

The  physical  picture  of  radiation  pressure  is  less  simple 
in  acoustics  than  in  electrodynamics,  where  we  have  to  do  with  the 
linear  set  of  Ilaxwell ' s  equations  and  with  transverse  waves.  In 
the  latter  case  the  radiation  pressure  is  determined  by  the  electric 
and  magnetic  field-strengths  and  the  formula  for  the  radiation 
pressure  is  valid  for  all  field  intensities. 

On  the  other  hand,  the  equations  describing  the  notion  of 
acoustic 'waves  are  in  general  not  linear.  To  avoid  mathematical 
difficulties,  the  equations  are  usually  "linearized"  by  using 
developments  of  the  non-linear  expressions  in  Taylor  series  and  re¬ 
taining  only  first-order  terms.  The  characteristic  quantity  of 
these  developments  is  the  particle-displacement.  The  results  of 
this  "linearized"  theory  are  therefore  first  approximations  which 
can  be  expected  to  be  valid  only  for  small  displacements.  Radiation 
pressure,  however,  is  connected  with  energy-densities,  which  are 
quadratic  terms,  containing  the  squares  of  displacements  or  veloci¬ 
ties.  Any  theory  dealing  with  radiation  pressure  must  therefore 
retain  at  least  all  second  order  terms  to  avoid  erroneous  results. 
Even  if  this  is  done,  the  results  are  still  only  valid  for  small. 
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amplitudes,  if  Taylor  expansions  are  used.  Simple  relations  between 
radiation  pressure  and  mean  total  energy  density  (as  for  instance 

*  rad  =  E)  are  found  only  for  small  amplitudes  of  displacement. 
Relations  of  this  kind  do  not  express  a  basic  law,  as  they  do  in 
electrodynamics,  where  they  are  independent  of  the  amplitudes  of 
the  fields. 

For  the  special  case  of  a  liquid  we  shall  show  that  for 
large  amplitudes  no  simple  relation  between  radiation  pressure  and 
mean  total  energy  density  exists.  Only  in  the  second-order  approx¬ 
imation  of  the  general  formula  do  we  find  proportionality  between 
these  quantities. 

Each  mass  particle  in  a  simple -harmonic  acoustic  compres- 
sional  wave  makes  sinusoidal  movements  around  the  point  it  would 
occupy  if  there  were  no  wave  notion.  The  customary,  and  mathematic¬ 
ally  the  most  convenient,  way  to  describe  the  wave  propagation 
makes  use  of  this  displacement  of  any  particle  from  its  original 
location.  All  other  physical  quantities,  as  velocities  or  pressures, 
are  related  to  the  moving  particles.  The  equation  of  motion  that 
relates  always  to  the  same  moving  particle  is  usually  called  the 
"Lagrangian  equation  of  motion”,  or  the  "equation  of  motion  in 
Lagrangian  coordinates". 

The  other  way  of  dealing  with  these  problems  uses  an  equa¬ 
tion  of  motion  that  relates  to  fixed  coordinates  in  space,  called 
"the  equation  of  notion  in  Eulorian  coordinates".  Here  we  are  not 
following  the  motion  of  a  certain  particle,  but  are  observing  what 
happens  at  a  fixed  point. 

A  careful  distinction  must  be  made  between  velocities. 


densities,  pressures,  and  other  quantities  related  to  the  instantaneous 
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position  of  a  moving  particle,  and  the  sane  quantities  when  related 
to  a  fixed  point  in  space.  In  dealing  with  the  mean  radiation 
pressure  upon  a  reflector,  for  instance,  we  shall  assume  the 
average  position  of  the  reflector  to  he  stationary  in  space.  But 
lie  may  also  speak  of  the  mean  pressure  observed  at  a  surface  which 

is  moving  together  with  the  particles  of  a  plane  wave* 

An  acoustic  coupressional  harmonic  t/ave  causes  two  impor¬ 
tant  effects:  1 )  It  changes  the  mean  hydrostatic  pressure  at  all 
points  affected  by  the  wave.  2)  It  causes  additional  mean  stresses 
in  the  medium  due  to  the  time  average  of  the  periodic  flow  of 
mechanical  momentum  of  the  compress ional  wave. 

If  a  plane  acoustic  wave  of  infinite  width  traverses  a 
medium  and  strikes  a  perfect  absorber,  the  mean  pressure  inside 
the  wave  region  undergoes  some  diminution  due  to  the  wave  notion. 

But  if  the  acoustic  beam  is  surrounded  by  undisturbed  regions  of  the 
medium,  as  is  ordinarily  the  case,  the  mean  pressure  becomes  equalized 
throughout  the  entire  fluid.  In  this  case  the  effect  1 )  mentioned 
above  does  not  cause  a  directed  force  upon  obstacles  placed  in  the 
path  of  the  wave.  The  basic  physical  cause  of  the  acoustic  radia¬ 
tion  pressure,  or  better  radiation  tensor,  is  in  this  case  the  time 
average  of  the  periodic  flow  of  mechanical  momentum  in  the  region 
affected  by  the  beam.  (See  Part  II,  Sec.  12). 

In  a  non-vlscous  fluid  affected  by  plane  acoustic  waves, 
there  is  no  net  flow  of  matter  due  to  the  wave  propagation.  (The 
unidirectional  "hydro dynani c  flow"  caused  by  acoustic  waves  is  an 
additional  physical  effect  due  to  the  viscosity  of  the  fluid  and  is 
not  to  be  taken  into  consideration  here,  see  also  Part  II,  Sec.  10). 

Considering  a  section  of  unit  area  of  a  plane  coupressional 
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wave  (Fig.  6)  which  we  imagine  as  fixed  in  space. 


-u(a,t)  f — S’-  u(a+da,t) 


x 


x=a 


x=a+t.a 


Fig.  6 


Particle  movement  through 
control-areas  fixed  in  space 

the  average  value  in  tine  of  the  mss  flow  is 


-  _  i  f  +  Tp 

r  -  vt  p 


(a,t)«  u  (a,t)  dt 


Tp  =  l/f  is  the  period, p  (a,t)  the  mss  density  at  any  moment  t  at 
the  fixed  coordinate  x  =  a  of  the  cross-section  under  consideration, 
and  u  (a,t)  the  instantaneous  velocity  of  the  particles  crossing  the 
section  a  at  the  time  t.  At  different  tines  different  particles, 
are  crossing  the  section  a,  so  that  u  andp  refer  only  to  the  particles 
that  cross  at  the  tine  t.  For  a  periodic  wave  notion  in  a  medium 
that  was  originally  at  rest  (u*=  0),  pu  vanishes  at  small  amplitudes. 

Ilext  we  consider  the  flow  (or  flux)  of  momentum  through 
the  section  a.  Arguing  in  the  same  way  as  in  the  preceding  case 
of  the  electromagnetic  flow  of  momentum,  the  momentum  mu  crossing 
unit  area  at  x  =  a  in  one  second  is  mu  =  (pm)  •  u  =  This 

quantity  is  obviously  always  positive.  This  means  that  the  plane 
wave  carries  periodically  momentum  -  whose  average  value  in  time 
does  not  vanish  -  through  fixed  sections  perpendicular  to  the  direc¬ 
tion  of  wave  propagation. 
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The  quantity  nu  =  (p  u)  •  u  -pu^  can  be  interpreted  in 
another  way,  leading  to  the  concept  of  the  "flux  of  momentum", 
pu  is  the  density  of  mechanical  momentum  and  u  is  the  velocity  with 
which  the  quantity  pu  crosses  a  section.  Thus  we  can  say  that  the 
"flux  of  momentum -density"  per  second  through  a  section  of  unit  area 

p 

equals  pu  •  u  =  p  u  ,  just  as  we  speak  of  a  "flux  of  mass-density", 
p *  u  per  second  over  a  unit  area. 

For  small  amplitudes  the  time  average  of  the  kinetic 
energy-density  l/2  p  u^  is  equal  to  the  time  average  of  the  poten¬ 
tial  energy.  In  this  case  l/2  p"u7  is  equal  to  E/2  (E  =  mean  total 
energy  density).  The  average  value  in  time  of  the  "flux  of  momentum 
density"  p u2  is  therefore  equivalent  to  the  mean  total  energy 
density  E. 

Considering  a  volume  element  between  two  unit  cross  sections 

at  a  and  a  +  da  (Fig.  6),  a  flux  of  momentum  p>(a,t)  *  u^(a,t)  enters 

the  section  a  and  a  flux  p{ a  +  da,t)  •  (a  +  da,t)  leaves  the 

^  (Cj  u*-} 

section  (a  +  da)  at  the  same  time.  The  mean  value  -  '  •  da 

O  a 

in  time  is  obviously  the  gain  in  momentum  of  the  volume  element 
between  a  and  a  +  da  per  second.  According  to  Newton 1 s  Law,  this 
gain  in  momentum  per  second  is  equivalent  to  a  force  exerted  upon 
the  volume  element.  The  volume  element  reacts  with  a  force 

+  pH  da  to  preserve  equilibrium. 


n  a 


lie  can  compare  this  picture  with  the  volume  element  tinder 
the  action  of  a  stress  in  the  x-direction  (Fig.  7). 
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a  a+cla 


Fig.  7 

Volume  element  under  stress 
in  the  x-direction 

--N,  m 

xx 

The  force  exerted  upon  the  volume  element  is  known  to  be  H - “  da. 

1  QQl 

2 

The  mean  flux  of  momentum  ^ou  can  therefore  be  considered  with  re¬ 
gard  to  its  mechanical  action  as  equivalent  to  a  mean  stress 

—  ~2 

=  -O u  ,  acting  on  any  cross  section  affected  by  a  compressional 
wave.  The  effective  mean  stress  component  T  in  the  fluid  due  to 
the  plane  wave  is  therefore 

^xx  eff.  ~  ^xx  P  u  ~  “  (p  + 

\  l 

p  is  the  mean  pressure  and  p  the  mean  flux  of  momentum  density 

l 

at  a  fixed  section. 

lie  have  attempted  to  present  a  simple  physical  picture  for 
the  case  of  a  plane  wave  that  is  in  accord  with  a  strict  theoretical 
treatment  of  the  motion  of  mechanical  waves.  In  the  general  case 
of  a  three-dimensional  notions  with  velocities  u^j  Uy,  uz?  it  can 
be  shown  that  the  resultant  stress  tensor  in  a  medium  in  Eulerian 
coordinates  becomes  (see  ref.  5S  pp.  241  and  290): 
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(1.15) 


The  components  of  the  flux  of  momentum  in  general  are  ( p 

and  are  physically  equivalent  to  the  components  of  a  stress-tensor. 

The  hydrostatic  pressure  of  the  fluid  in  absence  of  wave  motion  is 
denoted  by  pQ5  while  p  represents  the  change  in  pQ  due  to  the  wave 
motion  (excess  pressure).  The  total  effective  pressure  is 
Pt  =  Pn  +  p«  For  a  plane  wave  in  a  fluid  without  absorption ,  u  =  u  5 
Txx  ~  T-'/y  =  Tzz  ~  ~  Pt»  E<1*  (1.15)  reduces ,  in  time  average*  to 


-  (pt+  pu/  ) 

0 

0 

T  = 

0 

-  Pt 

0 

0 

0 

-  Pt 

(1.16) 


It  is  the  terra  (p  +  Cu  )  that  turns  out  to  be  responsible  for  the 

\  x 

radiation  pressure  exerted  by  a  plane  acoustic  wave.  According  to 

Eq.  (1.16)S  this  quantity  obviously  has  the  character  of  a  stress. 

Considering  a  material  surface  element  perpendicular  to 

the  direction  of  wave  propagation  x,  and  the  adjacent  volume 

element  of  a  fluid  (Fig.  3) ,  the  surface-element  must  exert  a  stress 

D 

_  *Pt|f 


Pt  + 


i!  ^  Po 

— H - ?■-  x 


xt 
Fig.  8 

Forces  exerted  upon  a  volume  element 
adjacent  to  a  device  D  for  measuring 
radiation  pressure. 
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-  (p.j.  +  pu,/^)  upon  the  right  side  of  the  volume  element  in  order 
to  maintain  equilibrium.  The  surface  therefore  undergoes  a 


pressure  p^  +  pu2  in  the  x-direction  on  its  left  side. 


The 


static  pressure  behind  the  surface  is  denoted  by  pQ.  The  result¬ 


ing  mean  pressure  exerted  on  the  surface  is  therefore s 


?rad  =  Pt  “  P0  +  PU*2  =  P  +  J°  \2 


(1.17) 


This  is  the  general  formula  for  the  radiation  pressure  of  a  plane 
wave  upon  a  plane  material  surface  perpendicular  to  the  direction 
of  wave  propagation.  If  the  wave  region  communicates  (for  ex¬ 
ample  by  a  small  hole  in  the  absorbing  surface)  with  the  medium 
behind  the  surface,  p^.  =  pQ  at  a  perfect  absorber  and  ~PUX2* 

For  small  amplitudes,  pu^  =  E,  and  therefore,  Prad  =  E.  (1.1  S) 

For  a  perfectly  reflecting  surface  the  physical  pic¬ 
ture  is  more  complicated,  as  we  shall  show  later  in  the  exact 
treatment^  here  p^  at  the  surface  differs  from  pQ  and  turns  out 
to  be  pQ  +  2E  for  small  amplitudes,  therefore  .  =  2E  for  such 
a  surface  perpendicular  to  the  wave  propagation. 

Thus  finally,  we  have  arrived  at  the  sane  relations 
as  we  found  to  hold  in  electrodynamics,  but  the  physical  back¬ 
ground  is  more  complicated  in  the  case  of  mechanical  waves. 

Other  cases  will  be  treated  mathematically  in  this  report.  Our 
purpose  here  is  to  give  a  preliminary  idea  of  the  general  causes 
leading  to  an  acoustic  radiation  pressure,  for  some  special  and 
simple  cases. 
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PART  II 


PLAI1E  COI IPRLSSIOIIRL  HAVES  AID  RADIATION  PRESSURE  IM  LIQUIDS 


1 •  General  Considerations  -  Compressibility 
As  a  whole,  this  report  deals  mainly  with  the  radiation 
pressure  of  plane  acoustic  waves  in  liquids.  There  are  two  special 
reasons  for  this. 


first,  because  this  case  is  of  practical  importance.  Al¬ 
though  a  plane  wave  cannot  strictly  be  realised  experimentally, 
still  if  the  width  of  the  acoustic  bean  is  large  in  comparison  with 
the  wavelength,  the  concept  of  a  plane  wave  gives  us  a  good  approxima¬ 
tion,  especially  in  the  case  of  the  high  frequencies  used  in  ultra¬ 
sonic  waves. 

Second,  because  in  the  acoustics  of  liquids  we  are  able 
to  make  use  of  the  concept  of  constant  compressibility.  This 
concept  introduces  a  simple  analytical  relation  between  the 
hydrodynamic  pressure  p  and  the  relative  change  in  volume  V/V0 
of  a  volume  element  having  the  original  volume  V0,  on  which  a 
pressure  Z\p  is  exerted.  The  compressibility  p  is  by  definition; 


A  t.  _i_ 

v0  Ap 

For  small  changes  in  volume  and  pressure,  ve  have  the  isothermal 
compres  sibility , 


Pa 


isoth. 


-L 

Vr 


(If) 

dp  ,.T 


(2.1) 


(2.2) 


and  for  the  adiabatic  compressibility 


(2.3) 


T9 


the  subscripts  T  and  <y  denoting  constant  temperature  and  entropy* 
For  acoustic  naves  one  uses  normally  the  value  s 

though  the  process  of  compression  and  rarefaction  of  the  liquid 
is  certainly  not  strictly  adiabatic,  owing  to  the  unavoidable 
dissipation  of  heat.  Still  the  exact  value  of  p  will  not  be 
very  much  affected  even  if  the  process  is  not  strictly  adiabatic. 
This  fact  can  be  seen  from  the  formula  giving  the  difference 


between  p  ,  and  £.  ,, 

^adiab*  ^ 1S0th. 

Padiab.  “  Pisoth. 


4. 


-! 

i 

Vr 


n>v\ 

\)pL\3p/_ 


Q)V 


(2.4) 


Numerically  it  turns  out  that  this  difference  is  rather  small 
for  liquids  forming  drops.  It  amounts  to  only  a  few  percent 
under  normal  conditions  of  temperature  and  pressure.  Thus  for 
water,  the  value  Padia-k,  =  46(10“^)  cm2/kg  force  at  t  =  20°G  is 
ordinarily  used  for  compress ional  waves. 

Over  a  large  range  of  pressure  the  compressibility  is 
not  constant.  At  t  =  20°C,  for  instance,  Pj_soth,  ^ias  ^le  follow¬ 
ing  mean  values  between  p  =  1  and  p  =  1  +  Ap. 


f> 

II 
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The  excess  pressure  in  a  plane  compressional  wave  is  a 
function  of  the  intensity  J  of  the  wave.  Its  maximum  value  is 
given  by 


P. 


max 


(2.5) 


where  c  is  the  velocity  of  sound. 
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For  water,  at  the  very  high  intensity  of  J  =  tOO  watt/cm^, 
p  amounts  approximately  to  17  atm.  Thus  over  the  range  usually 
encountered ,  the  compressibility  can  be  regarded  as  practically 
constant.  This  fact  greatly  simplifies  the  mathematical  treatment 
of  acoustic  waves  in  liquids,  and  enables  us  also  to  extend  the 
discussion  to  finite  amplitudes  in  liquids. 

2.  Strains  and  Stresses  in  Viscous  Liquids  for  Plane 
Congressional  Waves 

The  general  relations  between  stresses  and  strains  in  a 
viscous  medium  having  a  coefficient  of  viscosity  vj  and  a 
compressibility  £  are: 

rxx  =  -4-  +  |  h  <2  s  -  <v  -  s ) 

Py  =  |  1  (2  <y  -  S  -  Sc  )  (2-6) 


m  _ 

iZZ  ~ 

A 

P 

2 

■+  5  1 

(2  t  -  c  - 

Z  X 

cy  } 

li 

2i/|ct ys  Tyz  : 

2yjOyZ,  j.zx 

=  2  rl  £W 

A  = 

CX  +  cy  +  S 

c„  = 

C*^r  = 

22L  ,  = 

:  2A 

y 

97  z 

dz 

ii 

& 

o  1 

1 

“F 

(*£. 

3  y 

■^—7—)  etc, 

3  x 

£  = 

displacer 

rents  in  the  x 

;,  y,  z  directions 

A  dot  denotes  the  partial  time-derivative.  For  example, 

:  _  a.  c  =_a_  as  =  li. 

x  0  t  X  9  t  9  X  9  x3  ^ 
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For  the  sake  of  generality  we  include  here  viscous  forces s 
though  later  on  we  shall  neglect  viscosity.  The  concept  of  plane 
waves  means  that  at  all  points  of  any  plane  perpendicular  to  the 


direction  x  of  the  wave -propagation  the  state  of  motion  is  the 

same  and  therefore  —  =  -2-  =  09  also  the  displacement  Y}~  0  "  0. 

G »  o  2;  / 

Introducing  these  assumptions  in  Eq.  (2.6)  we  find: 


1  iJz  +  4v) 

P  3x  3  IdxBt 

i 

p  ^  3  /a*a t 


=  1  J2I  ’  2  y , 

P  2*  3  /^t 


(2.8) 


=  0  . 


The  pressure  p  is  defined  by: 

^xx  +  Tvy  +  Tzz 

“  P  =  - - 3 

The  quantity  T^  +  Tyy  +  Tzz  is  known  to  be  an  invariant  of  the 
tensor  T*j,  and  independent  of  the  orientation  of  the  axes  of  the 
coordinate  system.  This  is  the  reason  for  the  definition 
(2.9) s  which  is  also  in  agreement  with  our  earlier  definition 


A 

P 


(2.9) 


of  compressibility  in  Eq.  (2.1),  becau.se,  as  is  well  known 9  the 
quantity  /\,  =  Cxx  +  £yy  +  czz  ^-s  identical  with  the  relative 
change  in  volume  £\,  V/V0  of  a  volume  elenent  V0«  In  the  case 
of  a  nlane  wave  one  finds  from  Sqs.  (2.9)  and  (2.8): 


-  P 


=  iai 
P  3* 


(2.10) 


The  terms  containing  rj  cancel  out,  as  we  expect,  because  p  in 
(2.10)  gives  the  mean  pressure  in  a  volume  element  due  to  the 
wave  motion,  which  is  independent  of  viscous  forces  and  the 
same  in  all  directions.  In  a  plane  wave  a  volume  element  is, 
it  is  true,  stretched  and  compressed  only  in  the  direction  x 
of  the  wave  propagation,  because  the  displacements  have  only  one 
component  But  according  to  Eqs.  (2.1)  and  (2.10)  this  causes 
a  hydrodynamic  pressure  p,  which  is  a  scalar  end  therefore  in¬ 
dependent  of  direction.  The  relation  (2.10)  is  exact  by  defini¬ 
tion  for  media  with  constant  compressibility;  no  higher  order 
terms  in  £  need  be  considered. 

3.  The  Eulerian  Equations  of  Hot ion  for  Plane 
Comprossional  haves  in  Viscous  Liquids 

For  treating  the  radiation  pressure  in  a  plane  conpressional 
wave  we  must  know  the  stresses  in  the  medium.  As  the  radiation 
pressure  acts  on  the  surfaces  of  bodies  inserted  in  the  medium, 
whose  mean  position  in  space  can  be  regarded  as  fixed,  we  desire 
to  know  the  stresses  in  a  coordinate  system  referring  to  points 
fixed  in  space.  The  hydrodynamic  equations  applying  to  this 
case  are  the  so-called  "Eulerian1’  equations.  All  physical  quan¬ 
tities,  as  the  particle  velocity  u  or  the  pressure  p,  are  regarded 
as  functions  of  the  coordinates  x,  y,  z,  of  an  axial  system  fixed 
in  space,  and  of  the  time  t.  For  a  plane  conpressional  wave  with 
the  only  component  of  displacement  ^  in  the  x-direction,  the 
Eulerian  equation  of  notion  in  the  direction  x  of  wave -propagation 


is  known  to  be 
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where  p>  is  the  mass-density,  Du/Dt  the  total  derivative  with  re¬ 
spect  to  t  of  the  velocity  u  in  the  x-direction  of  a  volume -element 
having  the  density  p  .  For  the  plane  wave  here  considered, 
u  =  u  (x,t)  is  a  function  of  x  and  t  alone $  for  simplicity  we  use 
u  instead  of  u  .  Eq.  (2. 11 )  follows  immediately  from  consideration 
of  the  acceleration  which  a  volume- element  undergoes  under  the 

*■*!  rp 

action  of  the  force  in  the  x-direction.  No  accelerating 

forces  exist  in  the  y  and  z  directions,  because 


c_ 

$y 


dz 


=  0  and  Tik  =  0  for  i  f  lc 


From  (2o 11 )  we  have 


^  +  duto  _  ru,  „£u  = 

P[0 1  +  ax  at'  fit  +  2)x  dx 


(2.12) 


The  conservation  of  mass  of  a  volume  element  limited  by 
surfaces  fixed  in  space  (which  we  will  call  "control- surf aces")  re¬ 
quires  in  the  cne-dimensional  case  of  a  plane  wave 


-gig  +  ALpul  =  0 

at  +  e>x 

the  well-known  "equation  of  continuity". 

Combining  Eqs.  (2,12)  and  (2.13)  by  adding  u  4r t0 
(2.12),  and  using  ~^p^'  from  (2.13),  one  obtains  another  known  form 
of  Euler's  equation: 


(2.13) 


OX- 


cK  f>u)  +  <^(pu2)  _  ^Txx  _ 

s  t  ox  a  x 
Al gal  +  ALfpci  +  _ 


\ 

c-X 


(2.14) 


3t  ox 

This  is  the  main  equation  upon  which  we  base  our  further  considera¬ 
tion  of  the  stresses  in  a  plane  wave.  As  has  been  stated,  the 
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quantities  ^--'9  u,  and  are  functions  of  x  and  t,  where  x  is  a 
coordinate  fixed  in  spaces 


We  now  consider  a  volume  element  limited  by  two  planes  at 

x  and  x  +  dx  and  with  unit  area  perpendicu  lar  to  the  x-axis.  p  ,  u 

t 

and  T  are  tahen  with  respect  to  the  fixed  planes  x  and  x  +  dx 
respectively.  Wlien  particles  move  with  velocity  u  across  these 
planes,  u  (x,t)  belongs  to  different  particles  at  different  instants; 


that  is,  always  to  the  particular  particle  which  just  crosses  the 
stationary  plane  x  (or  x  +  dx)  at  the  time  to  p  is  the  density  of 
the  volume  element  between  the  "control-surfaces"  of  unit  area  at 
x  and  x  t  dx. 


We  can  interpret  the  terms  on  the  left  side  of  Eq.  (2.14) 
as  three  forces  acting  on  the  volume  element  between  the  control- 
surfaces  s  whose  sun  equals  zero  according  to  d'Alembert's  prin¬ 
ciple. 


The 


first  term 


d(iQu) 

3tT 


gives  the  change  in  mechanical  mo¬ 


mentum  |Du  at  x  with  respect  to  the  time  t*  This  can  be  inter¬ 
preted  as  equivalent  to  a  force  of  inertia  exerted  by  the  volume- 
element  5  as  is  usually  done  in  mechanics. 

The  term  -  -V-  -  means  the  mechanical  force  acting  on  the 
ox 

volume  element.  According  to  Eq.  (2.8)  it  is  given  by 


(2.15) 


The  third  term  is  due  to  the  gain  in  momentum  in 

C  x 

unit  time  which  the  volume  element  undergoes,  if  a  greater  amount  of 
momentum  pu  enters  the  area  at  x  than  leaves  the  area  at  x  +  dx 
per  second.  Paraphrasing  L.  Brillouin's  terminology  we  call  the 
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quantity  p  u  "the  flux  of  mechanical  momentum  density”  ("Flux  do 
quantite/ de  mouvenent";  see  ref.  5S  pp.  241 1  and  290).  It  plays  a 
dominant  part  in  dealing  with  radiation  pressure. 

The  incoming  flux  of  momentum  in  unit  time  at  x  is  obviously 
(^>u*u)x;  the  outgoing  flux  at  x  +  dx  in  the  same  time  is  ( ^  u'T-O^dx* 
The  difference  (  p  ^u^^  -  (|pu*u)x  is  given  by 

~s  /  2\ 

.  dx 
d  x 

This  quantity  is  a  contribution- to  the  momentum  of  the  element  in  just 

as  true  a  sense  as  the  term  ,  dx  and  can  be  regarded  there¬ 

at 

fore  in  the  same  way  as  a  force  acting  on  the  element. 

It  is  obviously  convenient  to  add  this  force  due  to  the  flux 


of  momentum  pv?-  at  x  to  the  stress -tensor  T, 


So  we  write  (2.114) 


thus : 


d>&u  d  ( m  _  n 

d't  “  3x  V 


u2)  =  0 


(2.16) 


The  flux  of  momentum  /O  u'"'  is  equivalent  to  an  additional 
stress -component s  acting  on  each  volume  element  of  a  medium,  the 
particles  of  which'  -are  in  motion.  This  stress  reduces  to  a  particu¬ 
larly  simple  term  in  the  case  of  plane  waves.  For  the  general  case 
of  a  movement  of  particles  in  different  directions  with  velocities 

u  s  u  s  u  in  the  xs  y  and  z-directions9  the  stress-tensor  T  is  to 
a  y  » 

be  completed  in  the  way  we  have  already  mentioned  in  (1,115), 

The  complete  dynamic  stress -tensor  in  Eulerian  coordinates 
for  a  plane  wave  is  given  by 
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T 

1xx  , 

2 

or 

T 

*y 

0 

T  = 

T 

xy 

T 

yy 

0 

0 

0 

T 

zz 

(2.17) 


with  the  values  T  as  given  in  Eq.  (2.8). 

Disregarding  viscous  forces,  we  find  by  introducing  the 
pressure  p  due  to  the  wave  notion,  according  to  Eq,  (2.10),  and 
adding  the  static  pressure  pQ  to  p: 


-  (p.  +  O  u2) 
t  1 

0 

0 

0 

-pt 

0 

0 

1 

0 

-pt 

(2.18) 


as  already  given  in  (1.1 6), 

Eq.  (2.17)  or  (2.18)  will  answer  all  questions  concerning 
the  radiation  pressure,  because  it  gives  the  stress-components  in  a 
plane  compressional  wave.  lie  see  that  the  medium  undergoes  a  non¬ 
isotropic  state  of  tension  due  to  the  unidirectional  flux  of  momentum 
in  the  x-direction. 

We  seek  the  solution  of  Eq.  (2.16)  for  the  case  of  pure 

o 

sinusoidal  waves  5  in  order  to  insert  the  values  of  p^.  and  in 

Eq.  (2.18). 

The  Eulerian  equation  (2.16)  or  (2.12)  in  fixed  coordinates 

3  u  c?  (  p  u2 ) 

is  nonlinear,  because  of  terns  like  u  ■~s"“  or  like  . 


This  difficulty,  which  complicates  the  solution,  can  -  at  least  in 
liquids  -  be  avoided  by  using  another  set  of  aquations  of  motion, 
which  are  usually  called  the  "Lagrangian  equations  of  motion".  These 
equations  relate  not  to  points  fixed  in  space,  but  rather  to  the  mov¬ 
ing  particles.  It  is  much  easier  to  find  the  correct  solution  by 
means  of  the  Lagrangian  equations?  for  this  reason  they  are  usually 
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employed  to  find  the  solution  for  plane  compressional  waves.  Never¬ 
theless,  for  expressing  the  radiation  pressure  exerted  on  obstacles 
whose  mean  position  in  space  can  be  regarded  as  fixed,  these  solutions 
must  be  transformed  from  the  system  related  to  moving  particles  into 
the  Eulerian  coordinate  system,.  This  can  easily  be  done,  as  will 
be  seena 

One  important  point,  which  must  hot  be  overlooked,  is  the 
circumstance  that  the  average  values  in  tine  of  both  terms  p  and 
p u  in  fixed  coordinates  are  second-order  terms.  The  radiation 
pressure  is,  as  has  been  seen,  proportional  to  the  energy-density, 
and  is  therefore  a  second-order  quantity^  this  explains  also  its 
relatively  small  numerical  values  in  comparison  with  first-order 
pressures,  even  at  the  small  amplitudes  ordinarily  used.  Uhereas 
the  first-order  pressure  has  maximal  values  up  to  some  kilograms  per 

p 

cm  ,  the  acoustic  radiation  pressure  only  reaches  values  of  the  order 
of  grams  or  dynes  of  force  per  cm^» 

Ir.  transforming  p  and  p  u^  from  the  Lagrangian  into  the 
Eulerian  system  all  second-order  terms  must  be  carefully  taken  into 
account.  Restriction  to  first-order  solutions  leads  to  erroneous 
results.  Instances  night  be  cited  in  certain  papers  dealing  with 
radiation  pressure. 

In  the  case  of  a  liquid  the  use  of  the  Lagrangian  equations 
gives  us  an  exact  s elution  for  finite  amplitudes,  which  can  be  trans¬ 
formed  into  the  Eulerian  coordinate  system.  In  the  case  of  gases, 
on  the  ether  hand,  even  the  Lagrangian  equation  becomes  nonlinear 
and  we  must  develop  the  Lagrangian  solution  in  series  at  least  up  to 
the  second-order  terms  in  order  to  find  the  radiation  pressure  for 
small  amplitudes. 
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4«  The  Lagrangian  Equation  of  I  lotion  for  Plane 

Congressional  Waves  in  Viscous  Liquids 

Let  £  be  the  displacement  of  a  particle  with  respect  to 

its  original  undisturbed  position  x.  If  there  is  a  wave  motion , 

particles  move  back  and  forth  through  x,  so  that  £  =  £  (xst). 

The  true  particle  coordinate  in  space  at  any  time  t  is  therefore 

x  +  £  (xst).  The  equation  of  notion  for  the  displacements  £  can 

9  12 

easily  be  established  and  is  found  to  be  5 


fc 


~N,2r  /  P}T* 

-  Z-3& 

dx 


at2 


This  equation  is  absolutely  exact,  with  nothing  neg¬ 
lected.  The  constant  is  the  original  undisturbed  density  in 
absence  of  wave  notion.  T*  (x5t)  is  the  stress  at  the  location 
of  the  moving  particle,  that  is,  at  x  +  ^  (xst).  Ue  use  the 
notation  T^.  for  stresses  associated  with  moving  volume  elements 
for  distinction  from  T  ,  in  bulerian  coordinates  related  to  points 
fixed  in  space.  The  particle  velocity  at  the  instantaneous  loca- 

^  r 

tion  x  +  £  (x3t)  is  denoted  similarly  by  u*(x,t)  =  g  t  * 

The  relation  between  the  quantities  T.^,  u,  and  in 
Eulerian  fixed  coordinates  and  the  same  quantities  T^.,  U" ,  and 
in  Lagrangian  coordinates  is  obviously  (see  also  ref.  16)  : 


T*  (x,t)  =  T  (x  +  ?,t) 
u*  (xst)  =  u  (x  +  £,t) 
p*  (x5t)  =  p  (x  +  £st) 


(2.19) 


(2.20) 
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or,  by  substituting  x  -  K  for  x, 

T*  (x  -  ?,t)  -  T  (x,t) 

u*  (x  -  £,t)  =  u  (x5t)  (2.21) 

P*u  -  K,t)  =p  (x,t) 


Eq.  (2.20)  or  (2.21 )  allows  us  to  find  T5  u,  and  p,  for  example, 
if  we  know  the  solution  for  T* ,  u*,  p#  in  Lagrangian  coordinates. 

For  the  special  case  of  liquids  with  constant  compressi¬ 
bility  (3  the  relation  between  T*  and  £  is  given  by  Eq.  (2.8),  and 
from  (2.19)  we  have  the  well-known  equation  of  motion  for  viscous 


liquids: 


1  ifl 

-\  o 
P  dx*- 


3  c)x2cH 


(2.22) 


This  equation  is  linear  and  can  therefore  be  solved  easily 
and  exactly,  whereas  the  corresponding  equation  for  the  same  case 
in  Eulerian  coordinates  is  nonlinear,  as  already  mentioned.  For 


gases  the  connection  between  T*  or  p*  and  the  displacement  in 

2QC 

general  is  nonlinear.  In  this  case,  the  solution  of  (2,1 j)  can 


only  be  expressed  in  the  form  of  a  series. 


5.  The  Exact  Solution  for  Plane  Compressional  haves 
in  lion-viscous  Liquids 

Since  we  shall  consider  radiation  pressure  here  without 
regarding  the  influence  of  viscous  absorption,  we  neglect  the  last 
term  on  the  right  side  of  Eq.  (2.22).  The  solution  for  plane 
waves  in  liquids  follows  at  once  from  this  equation.  For  pure 
sinusoidal  motion  of  a  particle  about  its  original  location  the 
solution  is  known  to  be 
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5  (x,t)  =  t 


Sin  (cot  -  kx)  +  Y  cos  (cat  +  kx  +  9)1  (2.23) 


cos 


with 


v  _  2n~  „  2iiT  _  co  —  oph>. 

k  -  -t~  -  -  -  —  ;  ca  -  27T1*, 

/V  c  c 


rr 


f°e 


ca 

k 


This  solution  corresponds  to  the  assumption,  commonly  made  with  re¬ 


spect  to  the  boundary  conditions,  that  the  wave  motion  is  generated 
by  a  piston-lilce  source  moving  harmonically  with  angular  frequency 
ca  around  its  average  position.  The  source  may  be  located  at  any 
x  or  at  infinity,  From  (2.23)  we  derive  the  following  quantities: 

u*  (x,t)  =  |  sin  (cot  -  lex)  +  Y_sin  (wt  +  lcx  +  9) )  (2.24) 

-p*(x,t)  =  =  j2  f~  =  ^{+0^  (mt  '  kx)  +  V-sin(rat  +  +  9)|(2-25) 

The  density  p#  follows  from 

pW«=poVQ  (2.26) 

V  is  the  original  volume  of  a  volume  element,  V*  the  volume  of  a 
o 

displaced  volume  element.  According  to  Eqa.  (2.1)  and  (2.10), 


V»  =  vo  +  V  =  V0  -  pvc  A  p« 

or  V*  =  V  (1  +  -i^)  (2.27) 

o  dx 

From  (2.26)  we  derive  the  density  p*  at  the  coordinate  of  the 
moving  particle: 

=  |  +  azU't)  (2,28) 

dx 

The  solution  given  by  Eqs.  (2.23,  2.24,  2,25)  are 
strictly  exact,  so  long  as  the  assumption  of  a  constant  compressi¬ 
bility  can  be  regarded  as  valid. 
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However,  there  is  one  important  restriction  which  limits 
these  solutions  to  a  definite  range  of  amplitudes  £  »  The  dif¬ 
ferential  equation  (2»22)  implies  the  supposition  that  £  (x,t)  and 
the  derivatives  u*  =  c)  ??/  dt  and  -pp*  =  d  are  unique  functions 

of  x  and  t.  At  a  coordinate  x.j  +  £  In  space,  for  instance,  only 
one  kind  of  particle  with  displacement  K,  (x^ )  and  velocity  u*  (x  ) 
at  a  certain  tine  t-is  assumed  to  exist  in  the  derivation  which 
leads  to  Eq.  (2.22),  At  the  moment  when  a  particle  A,  originally 
located  at  a  position  x^  behind  another  particle  B  originally  at  x^, 
undergoes  such  a  large  displacement  that  it  reaches  or  passes 

the  particle  B  at  x  +  ?(x  ),  we  would  find  two  different  particles 

B  B 

with,  in  general,  two  different  velocities  u*  (x.  )  or  u#  (x  )  at 

A  B 

the  sane  point  x^  +  £(x  )  =  x^  +  £(x^)  in  space  at  the  same  time. 

The  condition  to  he  imposed  upon  our  solution  (2.23)  for  £  for 
preventing  the  overtaking  of  one  partiole  by  another  is  found  to  be 


.21 


>-  i 


X 


This  can  be  understood  immediately  from  the  expression  for  in 

Eq.  (2,28).  If  two  different  particles  originally  located  at 

different  coordinates  x,  and  x„  cane  into  contact,  the  original 

A  B 

volume  between  the  planes  x  and  x  would  be  compressed  to  zero, 

A  B 

and  the  density  ^  would  become  infinite.  This  would  happen, 
as  may  be  seen  from  Eq.  (2.28),  if  the  denominator  vanished  or  if 
ciK/  c*x  -  -1.  Only  so  long  as  the  condition  above  is  valid,  does 
*  remain  finite;  from  this  fact,  together  with  Eq.  (2.23)  we  find 


f 


Oja 
<)  X 


-COS 


COS 


1  +sin  -  kx>  +  y -sin  +  kx  +  e>  >  >  -  1 


3a 


and  hence  the  condition  limiting  the  maximal  amplitude  of  any 
displacement  is 


<  1 


or 


A 

2? r 


(2.29) 


For  all  amplitudes  smaller  than  that  indicated  by 
(2.29)9  the  solution  (2,23)  is  exact  and  unique.  In  media  with 
constant  conpre s 3 ibility 5  conpressional  sinusoidal  waves  as  re¬ 
presented  by  Eq.  (2.23)  are  propagated  without  distortion.  This 
fact  has  in  principle  already  been  stated  by  Rayleigh,  though  he 
did  not  deal  specifically  with  real  liquids,  but  made  only  a 
theoretical  statement.  The  condition  (2.29)  is  more  a  theoretical 
than  a  practical  limitation  for  the  amplitudes  as  amplitudes  in 
the  neighbourhood  of  =  1  would  require  enormous  energies  that 
could  not  be  realized  experimentally.  (See  Sec.  8  b) 

By  means  of  the  transformation  formulas  (2.21)  we  are 
able  now  to  find  the  exact  solution  for  plane  compressional  waves 
in  liquids  in  Eulerian  coordinates.  Me  need  only  replace  x  in 
Eq.  (2.23)  etc.  by  x  -  £(x,t).  Choosing  the  sine  solution  for  £ 
in  (2.23)5  we  obtain  from  (2.24)  and  (2,25) 


(  -ip  — ,  1 

u  (x,t)  =  cn£0  jyos  i cot  -  lc  (x-0  j+  y  cos  j  cot  +  9  +  k  (x-0  j  j 

=  >  co.s  r  cot  -kx  +  k^Q  ^ sin  (cot  -  kx)  +y  sin  (cot  +  9  +  kx))j  ' 

+  y  cosjcot  +  9  +  lex  -  k^Q  ^sin(cot  -  kx)  (2.30) 


+  Ysin(cot  +  9  +  kx)) 


1 


and 


P  (x,t)  = 


k^o(  r  ,  ' 

"  1  cos  cot  -  kx  +  k£0  ’y  in  (cot  -  kx)  +  y  sin(cot  +  9  +  lc x)JJ 


,J3 


l 


-  y  cos  |  cot  +  9  +  lac  -  k^o  (yin  (cot  -  kx) 

-  •  +  y  sin(cot  +  9  +  kx)Jj 


(2.31) 
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These  expressions  (2.30)  and  (2.31 )  are  exact  solutions 

of  the  nonlinear  Eulerian  differential  equation  of  motion  (2. 14) 5 

for  amplitudes  £  <  >./ 2rr  according  to  (2.29)  and  for  the  special 
o 

case  of  pure  sinusoidal  waves  traveling  to  the  right  and  left 
given  by  the  Lagrangian  solution  (2.23)  for  the  displacement 
This  statement  can  be  verified  by  inserting  Eq.  (2,30)  and  the 
appropriate  expression  for  pi n  the  differential  equation  (2.14)« 
(See  Part  II,  Sec,  9  for  a  general  proof). 

The  independent  variables  x  and  t  enter  this  solution 
throughout  in  the  combination  cot  +  lex.  Thus  the  solution  in 
Eulerian  coordinates  also  represents  systems  of  waves  traveling  with 
the  phase-velocity  c  =  oo/k  to  the  right  and  left.  But  the  resul¬ 
tant  wave  notion  is  not  given  by  simple  superposition,  such  as  holds 
with  the  Lagrangian  coordinates.  Physically  this  means  that 
u(x,t)  and  p(x,t)  in  Eulerian  coordinates  are  not  represented  in 
our  case  by  superposition  of  pure  sinusoidal  traveling  waves*,  the 
wave  form  is  distorted  and  this  fact  is  responsible  for  the 
existance  of  higher  order  terns  in  ^0,  whose  average  values  in 
time  do  not  vanish. 

6.  Development  in  Series  up  to  the  Second  Order  of 

the  Exact  Solution  in  Eulerian  Coordinates 

For  snail  amplitudes  9  we  develop  the  solutions  (2.30) 

o 

and  (2.31)  in  Taylor  series,  retaining  second-order  terms  in 
For  u(x,t)  we  get  from  (2.30),  writing?!.  =  cot  +  lex: 

-r 

f  f  *v  **| 

u  =  o£o  I  cos  j  X  +  k9p  i^sin  X>  +  ^  sin(  Z'+  +  ®)/  I 

+  ^cos  I  +  9  -  lc90  sin  X  +  "'|,'sin(  t?+  +  9)^5 
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u  ~  °^o  |^cos^_  "  k^o  (ssin  +  y  sin(  '£+  +  Q)J  *  sin  T 
+  ycos(£j.  +  9)  +  y  k£0  ^  sin  77  +  y  sin  (7^.  9) y 


sin(  '£'  +  9) 


(2.32) 


As  we  shall  be  chiefly  interested  in  average  time  values 
of  u  and  p5  ue  compute  the  mean  value  u  of  u  by 


u(x,t)  = 


t,  +  T 
1  I  1  P 


u(x,t)  dt  = 


„cjot„  +  oiT 

1  1  P 


u(xscnt)  dcot 


.t  +  2rr 

_L  (  1 

u(xst)  =  2ir  J  u(x,oot)  dcot 


(2.33) 


From  (2.32)  and  (2.33)  ue  find 

2.  2 

u  =  —  (l  +  ycos(2  lex  +  9)  )  +  — (y  +  cos (2  lex  +  9)  ) 


o <k.Z  ,  '  2^ 


(2.34) 


The  result  (2.34)  shows  that  u(xst)  does  not  depend 

upon  x  even  if  there  is  a  reflected  wave,  the  latter  being 

characterized  by  the  reflection  coefficient  y .  In  Lagrangian 

coordinates  the  tine  average  u*  of  the  velocity  u*  equals  zero. 

In  the  sane  way  as  in  deriving  u5  we  find  from  Eq. 

(2,31)  for  the  average  tine  value  of  p, 

2  2 

_  k  l  2 

p  - - —  (1+2  cos (2  kx  +  9)  +y  ) 


(2.35) 
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The  mean  pressure  p  in  Eulerian  coordinates  evidently 
varies  sinusoidally  with  x  when  reflected  waves  are  present. 

Another  mean  value  that  we  shall  need  is  the  time-average 
value  of  the  quantity  ^u^  in  Eulerian  coordinates.  If  we  trans¬ 
form  p*  into  p)  and  develop  p  in  powers  of  tQ  ,  we  obtain 

fo  (1 +ai?o  +  a25o2  +  —) 

The  development  of  u  may  be  written 

u  =  b  £  +  bj^’2  +  ...  =■  l  (b  +  b  l  +  ...)  (2.36) 

10  2o  ol  2  o 

as  seen  from  (2.32).  So  we  find 

f  “2  (1  +  alh  +  a/o2)  '  (b1  +  V-o>2  ?o2+  — 

or,  disregarding  terms  of  higher  than  the  second  order, 

pu2  =  b2  £  2  +  ....  =  c*u*2  (2.37) 

Therefore  ire  do  not  need  the  development  of  p  in  powers  of  as 
only  the  undisturbed  density  p>o  enters  (2.37).  This  relation 
shows  also  that  the  value  of  the  kinetic  energy  density  at  small 
amplitudes  is  the  sane  in  both  Eulerian  and  Lagrangian  coordinates. 

We  have  from  Eqs.  (2.36),  (2.37)  and  (2.32), 


pu2  _  po£o2  .  =  p0^02  *  Cco  cos  ?r_  +  co  ycos(  +  ©)  )2 

or 


(1  +2|  cos  (2  kx  +  9)  +  Y2)  = 


3 


f 


*2 


u 


(2.38) 


Comparing  p  and  pu2  in  (2.35)  and  (2.38),  we  find  p  = 
This  relation  holds  in  liquids  at  small  amplitudes,  but 
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not  in  gases.  (See  Part  III,  Sec.  4). 

It  is  conformable  to  our  purpose  to  introduce  the 
mean  total  energy  density  E  in  (2.35)  and  (2.33).  For  a  pure 
sinusoidal  plane  compressional  wave  traveling  to  infinity  in 


one  direction,  we  get  the  mean  kinetic  energy  density  -p- 


u 


2 


directly  from  (2.33),  if  we  put  ^  =  0.  Therefore  for  small 
amplitudes 


2  ^2 
CO  o, 


E,  .  =  — t— 

Kin  4 


o^o 


E"- 


kin 


The  identity  of  Fj^n  in  Eulerian  coordinates  at  small  amplitudes 


with  E*^^  =  in  Lagrangian  coordinates  has  already  been 

seen  in  Eq.  (2.37). 

The  -potential  energy  density  due  to  the  elastic 
compression  of  the  medium  is  given  by  dE^^  =  -  p  •  dV,  that  is 
the  work  done  upon  unit  volume  VQ  =  1i  under  the  action  of  a 
pressure  p.  From  Eq.  (2.1)  we  find,  with  VQ  =  1 ,  dV  =  -  p  dp 

and  therefore  dE^^  =  p  pdp.  By  integration  we  have 

E*Pot  =  5-P  P*2 

Inserting  p#  from  (2.10)  we  obtain 

E*  -1  ii2 

L  pot  “  2p 

Since  according  to  Eq.  (2.23),  is  a  function  of  the 
argument  cot  +  lac  alone, 


II  =  t  iE  II  =  m  1  M  =  x  ^ 

dx  co  $t  c  c 


So  we  have,  remembering  that  d  -  i/p  n. 


(2.39) 
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E*. 


1  u*' 


pot 


2£ 


r  c 


2  = 


1  2 

2  pou" 


(2.40) 


The  difference  between  n  and  E#  .j.  in  Lagrangian  coordinates 
is  therefore  E#j_in  -  E*^^  =  (jO*  -  f>0)  a*2/2  in  liquids.  It 
vanishes  at  small  amplitudes s  that  is  in  the  second -order  approxi¬ 
mation. 


At  small  amplitudes  the  average  time  value  for  the 
potential  energy  density  in  pure  sinusoidal  waves  is  given 
therefore  by  the  same  expression  as  for  the  mean  kinetic  energy 
density  in  (2.39); 

.,,2  n  r  2 

E*  .  =  —  =  E*  .  (2.41) 

pot  ‘  4  xxn 


The  total  mean  energy  density  of  a  plane  unidirectional  wave  at 
small  amplitudes  is 


E#  =  E*,  .  +  E*  . 

kin  pot 


2  r  2 

{3qZo 


E 


(2.42) 


For  small  amplitudes  there  is  therefore  no  distinction  between 
E*  and  E. 

Introducing  Eq.  (2.42)  in  (2.38)  and  (2.35)  and  re¬ 
membering  that  k2/P  =  co2/pc2  =  p0(n25  we  find 

p  =  -  E  (1  +2  V'cos(2  kx  +  9)  +  ^  2)  (2.43) 


+  E  (l  +  2  y  cos (2  kx  +  9)  + 


(2.44) 


I-fe  :have  now-  derived  the  expressions  in  Eulerian  coordinates  needed 


to  establish  the  stress  tensor, 


30- 

O 


7.  Stress- tensor  and  Radiation  Pressure  for  Small 
Amplitudes  in  Liquids 

By  inserting  Eqs.  (2.43)  and  (2.44)  In  (2.18) ,  we  obtain 
the  time  average  of  the  stress  tensor  T  in  a  liquid  traversed  by 
plane  compressions!  waves  with  small  amplitudes  (k^o«1i)3  dis¬ 
regarding  viscous  absorption.  As  the  liquid  is  under  an  additional 
hydrostatic  pressure  p  ,  we  must  add  this  pressure  to  the  hydro- 
dynamic  pressure  p  in  (2.43)  due  to  the  wave  motion  alone.  So 
we  have  for  the  mean  stress  tensor  T 


-p°  ■  ...  ° 

0 

T  = 

0  -(pQ-  E[1  +  2  ycos(2kx  +  9)  +  y21) 

0 

0  0  -(pQ-  Ejj  +  2  v 

/  cos(2kx  +  9)  +y2]) 

(2.45) 


Owing  to  the  fact  that  the  two  quantities  p  in 
(2.43)  and  pu2  in  (2.44)  are  numerically  equal  but  of  opposite 
signss  the  sum  p  +  pu2  equals  zero.  Thus  in  (2.45)  is  equal 
to  the  static  pressure  p0  and  is  not  changed  by  the  wave  motion. 
The  existence  of  a.  concessional  wave  evidently  changes  only  the 
components  T,~.  and  T„^  perpendicular  to  the  direction  of  wave- 

<J  t  f  J  a 

propagation.  From  Eq.  (2.43)  it  is  seen  that  the  presence  of  a 


compressional  wave  in  a  liquid  diminishes  the  mean  static  pressure 
p0  by  the  amount  E  (l  +  2  cos (2kx  +  9)  +  ^  2),  In  the  direc¬ 
tion  x  of  wave  propagation  this  diminution  is  exactly  compensated 
by  the  stress  -pi  u2  due  to  the  f lux  of  momentum  in  the  same  direction. 

Now  let  us  insert  -  an  infinite  plane  material  surface 


perpendicular  to  x.  The  compressional  wave  traveling  to  the 
right  undergoes  reflection  at  this  surface,  characterized  by  the 


coefficient  y  of  tlie  reflected  amplitude  and  the  phase  angle  9 
of  the  reflected  nave.  For  computing  the  force  acting  upon 
the  reflector,  we  assume  that  behind  the  reflector  is  the  same 
liquid  under  the  static  pressure  p  •  The  resultant  force  per 
square  centimeter  acting  upon  the  reflector,  which  is  the  radra1 
tion  pressure,  equals  the  difference  "between  the  pressures  on 
the  left  and  right  sides  of  the  reflector.  (Fig.  9) 


Fig.  9 

Plane  reflector  undergoing  &.  pressure  pqon  the 
front  side  and  a  pressure  p0  on  the  rear 

The  pressure  exerted  at  the  boundary  between  the  reflector  an^  tne 
irradiated  medium  is  given  by 


Xn  =  T^  cos (n,x)  +  T^y  cos(n,y)  +  Txz  cos(n,z)  (2. 

Xn  denotes  the  component  of  pressure  in  the  x-direction 
and  n  the  direction  of  the  inner  normal  to  the  surface.  In  our 
case  (Fig.  9)  only  cos(n,:c)  =  -  1  is  different  from  zero,  further¬ 
more  T  r  =  T  =0  according  to  (2.45).  There  remains  only  a 
pressure  X„  =  -  T^  =  pQ  perpendicular  to  the  reflector  on  its 
left  side.  As  the  same  pressure  pQ  exists  at  the  right  side  of 
the  reflector,  no  resultant  force  is  exerted  upon  the,  reflecting 
surface.  A  plane  compresslonal  wave  extending  to  infinity,  in 
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would  exert  no  radiation  pressure  upon  any  plane  infinite  reflector 
perpendicular  to  the  direction  of  the  wave  propagation. 

This  case  of  a  plane  wave  extending  to  infinity  is  merely 
theoretical,  as  it  can  hardly,  even  approximately,  be  realized 
experimentally.  The  acoustic  beam  always  has  a  finite  cross- 
section  and  is  usually  surrounded  by  a  part  of  the  same  liquid  which 
is  not  affected  by  the  wave  motion.  The  mean  pressure  in  the  sur¬ 
rounding  liquid  is  p^j  the  pressure  inside  the  beam  is  changed  by 
the  wave  motion  to  pQ  +  p.  As  p  is  negative  according  to  Eq.  (2.43 )s 
the  mean  pressure  inside  the  beam  is  lox-rer  than  in  the  surrounding 
medium. 

The  mean  pressure  tends  to  be  equalized  over  the  whole 
liquid.  This  means  that  the  surrounding  liquid,  where  the  pressure 
pQ  is  assumed  to  be  maintained  constant,  compresses  the  beam-region 
until  the  mean  pressure  in  the  beam  is  the  same  as  that  around  it. 

This  effect  leads  to  a  radiation  pressure  upon  a  finite  material 
obstacle  placed  in  the  way  of  the  acoustic  beam. 

The  simplest  case,  and  that  which  has  received  most 
attention  in  the  literature,  is  the  radiation  pressure  upon  a  per¬ 
fectly  absorbing  surface.  There  no  reflected  wave  exists  and 
=  0,  The  theoretical  concept  often  used  to  represent  a  per¬ 
fect  absorber  regards  the  absorbing  wall  as  free  to  move  in  a 
direction  normal  to  its  plane  and  following  the  movement  of  the 
particles  immediately  adjacent  to  the  wall.  The  pressure  observed 
at  the  wall  is  identical  then  with  the  pressure  p*  associated  with 
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the  moving  particles,  which  \ie  knew  to  be  purely  sinusoidal  fren 
Eq,  (  2*  25),  whether  the  beam  is  finite  or  infinitely  extended. 
The  mean  work  done  by  p*  per  second  and  square  centimeter  at  the 
wall  equals  the  average  value  in  time 


1  rp  H 

-  t>« 


P*  &>  =  m 


T 


Pw  dt 


dt, 


which  by  use  of  Eqs.  (2.24;  and  (2.25)  with  Q,  and  from 
(2.42),  is  easily  computed  to  be  equal  to  E  •  c.  Such  a  moving 
wall  would  absorb  indeed  all  the  energy  reaching  it,  and  there 


would  be  no  reflection.  The  mean  force 


1  f: 

t„  *L 


P  o 


P"  dt  exerted 


upon  such  a  moving  wall  would  be  aero,  as  p*  is  purely  sinusoidal. 
The  freely  moving  and  absorbing  wall  would  be  indeed  a  perfect 
absorber,  but  it  represents  more  a  theoretical  fiction  than  a 
feasible  experimental  device. 

In  measuring  radiation  pressure  we  do  not  follow  the 
movement  of  the  particles  of  the  surface  struck  by  the  acoustic 
beam.  All  particles  on  the  surface  move  periodically  around 
their  original  positions,  so  that  the  center  of  mass  of  the  sur¬ 
face,  when  averaged  over  a  whole  period,  can  be  regarded  as  fixed 
in  space.  Usually  the  surface  subjected  to  radiation  pressure 
is  also  connected  to  a  measuring  device  of  considerable  inertia, 
■unaffected  by  the  rapid  motion  of  the  particles  in  the  wave. 

This  is  the  reason  why  we  introduce  the  Eulerian  system  of  coordi¬ 
nates  fixed  in  space  for  computing  the  time-average  of  the 
radiation  pressure  related  to  a  device  which  is  assumed  to  be 
fixed  in  space  rather  than  to  follow  the  motion  of  particles. 

(See  also  Sec.  12). 
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As  in  optics s  the  most  practicable  approach  to  a  perfect 
absorber  is  the  "hohlraum"  or  radiation  trap.  An  acoustic  hohlraum 
consists  of  a  cavity  with  acoustically  insulating  walls,  filled  with 
an  absorbing  medium,  and  provided  with  a  small  aperture  through 
which  the  acoustic  beam  is  admitted.  Such  a  device,  in  the  form 
of  a  cylindrical  tube  closed  at  one  end,  has  been  used  for  measuring 
acoustic  intensities  in  water.*  For  frequencies  in  the  megacycle 
range  the  absorption  of  energy  is  practically  complete  in  a  tube 
that  is  not  excessively  long.  The  plane  of  the  aperture  therefore 
serves  as  a  totally  absorbing  surface. 

Ue  now  consider  first  the  radiation  pressure  exerted  upon 
a  perfectly  absorbing  surface,  and  then  the  somewhat  more  complicated 
case  of  a  reflecting  surface. 

a.  Radiation  pressure  upon  a  perfectly  absorbing  surface. 


of  T: 


* 


(=  0 


2»45 ) 

we  find,  for 

0 

0 

II 

the  average-time  value 

“Pc 

0 

0 

0 

-(p0-E) 

0 

(2.47) 

0 

0 

-(p0-E) 

As  has  already  been  pointed  out,  the  mean  pressure  in  the 
acoustic  bean  is  lowered  by  the  wave  motion  by  an  amount  equal  to 
the  mean  energy  density  E.***  If  the  beam  of  infinite  cross-section 


*  Contract  I'l6  0KR-262,  Task  Order  Is  Report  !!o.  1,  Dec.  20,  1947, 
p.  36;  Report  Ho.  3S  Feb.  21,  1949s  pp.  8-11. 

**  See  also  ref.  5,  pp.  302-304 

****  An  instructive  physical  outline  of  this  effect  is  given  in  a 
paper  ‘03*  Hertz  and  Mende°.  (See  also  ref.  T') 
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does  not  communicate  with  the  region  of  the  liquid  unaffected  by 
the  wave  motion,  the  radiation  pressure  upon  an  infinite  wall, 
one  side  of  which  is  struck  by  the  wave  and  whose  other  side  is 
adjacent  to  the  undisturbed  liquid,  is  zero  because  both  sides 
of  the  wall  are  subject  to  the  same  force  pQ  per  square  centimeter 
in  the  x-direction. 

The  change  in  the  time-average  value  of  the  density  * 
due  to  the  wave  motion  is  found  from  Sq.  (2.28),  For  small 
amplitudes,  regarding  p*  as  a  small  quantity,  we  have  = 
po  (Ti  +  pp#  +  p2p*2  +  ...).  After  transformation  into  Eulerian 
coordinates,  \re  obtain  the  mean  change  in  density  at  small 


amplitudes 


.  Af  =  f  -  f0  ^  foP  (p  +  Pp  ) 


(2.48a) 


As  the  pressure  inside  the  beam  is  lowered  to  pQ  -  E,  p  equals  -E. 

pp^  is  found  immediately  from  Eqs.  (2.31)  and  (2.42)  to  be  +  E 

for  ^=0  and  small  amplitudes.  Therefore  we  obtain  from  (2.48a) 

7s,  p=  0,  which  means  that  the  wave  motion  in  the  infinitely  ex- 
\ 

tended  beam  does  not  change  the  mean  density  in  the  second-order 
approximation,  though  it  lowers  the  mean  pressure  by  the  amount  E. 

If  the  beam  is  of  finite  cross-section  and  surrounded 
by  liquid  with  hydrostatic  pressure  p0,  the  pressure  inside  the 
beam  being  pQ  -  E,  the  liquid  in  the  beam  undergoes  a  compression, 
until  the  mean  pressure  inside  reaches  the  same  value  pQ  as  in  the 
surrounding  region.  The  change  in  density  /\  p!  due  to  this 
compression  is  found  from  £\  n'/ ^)0  =  -  A  V'A0  =  P  A  P  or 


Ap'  =  p0P  Ap  = 


_  ^  p  _  E 


(2.48b) 
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Thus  the  wave  notion  causes  a  slight  increase  in  density  A^P* 
inside  a  bean  of  finite  width,  though  the  pressure  becomes  the 
same  in  the  bean  as  in  the  surrounding  undisturbed  medium  with 
the  unchanged  density* 

This  fact  can  perhaps  best  be  understood  by  the  follow¬ 
ing  reasoning:  The  acoustic  wave  motion  changes  the  pressure  by 
a  second-order  terra  (-E),  whereas  the  change  in  density  is  small 
of  higher  order.  The  compression  of  the  beam  of  finite  width  by 
the  surrounding  medium  changes  both  pressure  and  density  by  second 
order  terms  (E  and  E/c2  respectively).  In  the  resulting  pressure 
the  tiro  changes  cancel  each  other  to  zero  so  that  the  final  pressure 
becomes  pQ  again;  in  the  final  density  however  the  sum  of  the  tiro 


changes  results  in  the  second  order  term  /\p'  =  E/c2. 

I 

The  resultant  stress  tensor  becomes  now  (see  also 
ref.  5S  P«  302) 


“(P0  +  E) 

0 

0 

T  = 

y=e 

0 

"Pc 

0 

0 

0 

-p0 

(2.49) 


According  to  (2.40)  a  resultant  force  E  per  square  centimeter  is  now- 
exerted  upon  a  fixed,  perfectly  absorbing  wall.  The  radiation 
pressure  upon  such  a  surface  therefore  equals  the  mean  energy 
density,  as  long  as  the  amplitude  is  small: 


^rad)^-_  0  E 


(2.50) 


*  See  also  the  paper  by  G.  Richter^,  who,  following  a  different 
path,  arrives  at  the  same  result  concerning  the  mean  density 
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The  radiation  pressure  in  this  case  is  due  to  the  mean  periodic 
flux  of  mechanical  momentum  due  to  the  wave  motion,  which  intro¬ 
duces  the  term  ^Ou^  =  E  in  T^.  It  is  independent  of  the  direc¬ 
tion  of  propagation  of  the  wave  (as  is  expected  because  it  is  a 
tensor-component  and  not  a  vector) ;  this  means  that  Pra(j  has  the 
same  direction,  whether  an  incident  beam  strikes  the  absorber  or 
a  reflected  bean  leaves  it.  The  radiation  pressure  is  not  caused 
by  the  increase  in  density  ^  p'  inside  the  beam,  due  to  compres- 

g  1 

sion  of  the  surrounding  liquid,  as  has  often  been  asserted.  9 

In  the  plane  wave  of  infinite  width,  not  in  communication 
with  undisturbed  regions  of  the  liquid,  the  term  T^,  responsible 
for  the  radiation  pressure,  was,  according  to  (2.18),  found  to  be 
pQ  -  E  +  pu2.  The  two  terras  -E  and  pu2  =  +  E  cancel  exactly, 
which  means  that  the  diminution  in  mean  pressure  is  just  compen¬ 
sated  by  the  flux  of  momentum.  In  the  x-direction  no  additional 
stress  acting  on  a  volume  element  of  the  liquid  would  be  found, 
if  this  were  all  that  happened.  By  virtue  of  the  interaction 
with  the  surrounding  liquid  the  mean  pressure  pQ  -  E  is  brought 

up  to  the  value  p0  of  the  undisturbed  liquid.  The  term  T^  = 

—  2  2 
pQ  -  E  +  pu  is  changed  into  pn  +  pu  .  how  an  additional  stress 


pu*-. 


u  acts  throughout  the  beam  in  the  x-direction  and  leads  to  a 

1 


e 

radiation  pressure  E  =  pu^  upon  a  perfectly  absorbing  wall  per¬ 
pendicular  to  the  direction  x  of  wave  propagation. 

b.  Radiation  pressure  at  a  reflecting  surface. 
disregarding  viscosity 

In  the  case  of  a  reflected  wave  the  mean  pressure  over 
the  beam  is  given  by  Eq.  (2,43)  us  pn  -  E(l  +  2  ^cos(2kx  +  0)  +Y  ). 


a 
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It  varies  sinusoidally  along  the  x-axis.  If  the  beam  were  in¬ 
finites  no  radiation  pressure  would  result  on  an  infinite  sur¬ 
face,  as  has  already  been  pointed  out.  The  change  in  the  mean 
pressure  is  exactly  compensated  at  any  x  by  the  flux  of  momentum 
given  by  Eq.  (2.44).  The  resultant  stress  equals  pQ  as  in 
the  case  where  y  =  0. 

If  the  bean  is  of  finite  section  traversing  a  liquid 

with  hydrostatic  pressure  pQ  the  surrounding  liquid  still  has  the 

tendency  to  raise  the  pressure  inside  the  beam  to  pQ»  But  a 

uniform  pressure  ]^=  p  along  the  beam  is  not  compatible  with  the 

differential  equation  (2.14)  of  our  problem,  except  for  the  special 

case  in  which  y  =  0  (Sec.  7a).  This  can  also  be  seen  from  the 

solution  for  p  in  Eq.  (2.43).  When  a  reflected  wave  exists, 

p  necessarily  varies  periodically  in  space  along  the  x-axis.  It 
"0 

is  reasonable  to  conclude  therefore  that  by  action  of  the  surround¬ 
ing  medium  the  average  value  in  space  of  the  total  pressure 

p  =  p  +  P  is  brought  to  p  ;  inside  the  beam  o  varies  periodically 
to  o  *  t 

along  the  x-axis  around  the  new  value  pQ» 

At  points  remote  from  the  beam  the  pressure  in  the  un¬ 
disturbed  medium  is  assumed  to  be  homogeneous  and  equal  to  pQ. 

Within  an  "edge-region”  of  the  beam  a  transition  takes  place  from 

the  periodically  varying  pressure  p  ,  which  is  postulated  to  exist 

t 

in  the  beam,  to  the  homogeneous  static  pressure  pQ  in  the  surround¬ 
ing  medium.  In  this  transient  zone  a  more  complicated  (rotational) 
motion  of  particles  will  occur.  A  closer  theoretical  investigation 
of  this  effect  is  beyond  our  scope;  but  we  may  reason  that  the 
actual  transient  zone  will  be  small  if  the  wavelength  inside  the 
beam  is  a  small  fraction  of  the  width  and  length  of  the  beam. 
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Our  conclusion  that  the  space- average  in  p^.  inside  the 
beam  will  become  identical  with  the  static  pressure  pQ  outside  the 
beam  (strictly  at  y  =  z  =  <*  ) s  can  be  based  on  the  following 
theoretical  consideration:  Taking  the  time  average  of  the  differen¬ 
tial  equation  (2.14)s  and.  assuming  a  periodic  solution  in  time. 


we  get 


or 


ik  (p +  <^2>  =  (p +  fu2)  =  0 

p  +  pvfi  =  Constant 
\ 


(2.51.) 


This  constant  is  obviously  independent  of  t  and  x.  But  at  small 

amplitudes  it  is  also  independent  of  the  amplitude  of  the  wave 
motion,  or  of  the  quantity  k£Q,  since  the  sum  (p  +  pu^}is  independent 
of  k£o  as  seen  by  adding  Equations  (2.43)  and  (2.44)*.  In  the  case 
of  a  plane  wave  of  infinite  width  we  find,  the  constant  to  be  zero. 

For  a  finite  bean  within  an  undisturbed  medium,  however,  the  constant 
becomes  E  (1  +  y2).  Indeed,  if  we  assume  that  the  amplitude  is 
gradually  decreased,  that  is  kgQ  0  or  u->0,  the  total  pressure 
Pt  averaged  in  time  and  space  along  the  x-axis  must  tend  towards  pQ, 
the  static  pressure  in  the  surrounding  medium.  Therefore  for  u  -9  0 

A  +  X 

i.l  =  p 

3. 


o  +  X 


a  +  A 

p  dx  ->p 


or 


Lp4-K) 


-II  \ 

j  *  ■'a 


„a  +  A 


p  dx 


=  0 


u-^0 


*  See  also  Part  II,  Sec.  9 
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From  Eq.  (2.43)  we  find  that  we  have  to  add  the  constant  E  (l  +  \j  ) 
to  p  in  order  to  make  its  mean  value  in  space  vanish.  (Any  addi¬ 
tive  constant  in  p  is  compatible  with  the  solution  of  Eq.  (2.14)). 
Thus  under  the  action  of  the  surrounding  medium  the  pressure  p 
becomes 


-  E  (1  +  2  y  cos(2kx  +  9)  +  ^2)  +  E  (1  +  y  *) 


or 


P 


¥ 

2E  ycos(2kx  +  S) 


y 


u-^  0 


The  Constant  =  p  +|Or  now  reaches  the  value  E  (1  +  \^),  as  seen 

from  Eq.  (2.44) 5  in  agreement  with  our  statement  above. 

Since  the  constant  is  independent  of  the  amplitude  £0, 

the  foregoing  conclusion  holds  not  only  for  u  ->  0,  but  for  any  (small) 

.  fa  +  A_ 

value  of  u.  The  mean  total  pressure  ^  j  p^  dx  is  therefore 


a 


equal  to  pQ;  that  is  the  mean  value  of  p  in  space  is  zero,  in  the 
case  under  consideration. 

The  resultant  stress-tensor  is  obtained  therefore  from 
Eq.  (2.45)  by  increasing  the  mean  pressure  in  space  by  the  amount 
E  (1  tNj'2)  : 


T= 


•(po+  E(1  +\n) 


^  ))  0  0 

0  -(p0-2yl  cos(2kx  +  0))  0 

o 

0  0  -(p-  2^  E  cos(2kx  +  9)) 

(2.52) 

The  stress  tensor,  averaged  in  time  and  space,  becomes 


— 

-(p0  +  E(1  +  ^2)) 

0 

0 

T  = 

0 

"Po 

0 

0 

0 

-Po 

(2. 53) 


This  leads,  using  the  sane  considerations  as  above,  to  a  radiation 
pressure  upon  a  reflecting  surface,  whose  reflection  coefficient 
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is  y,  of  the  amount 


P 

rad 


(2.54) 


For  a  perfectly  reflecting  surface 
result 


Ti)  we  have  the  well  known 


'  ^ 
rad J 

y  =  i 

u 


-  Z  E 


(2.55) 


In  the  general  case  of  a  reflecting  surface  characterized 
by  the  amplitude  reflection  coefficient  y  and  the  phase  angle  0 
of  the  reflected  wave ,  the  radiation  pressure  is  independent  of 
this  angle  9.  Here  too  only  energy  densities  are  involved,  gs 
Y  ^  E  represents  the  mean  energy  density  of  the  reflected  wave 
component.  The  resulting  equation  (2.54)  can  therefore  be  described 
by  saying  that  the  radiation  pressure  is  composed  of  two  parts; 

One  part  is  due  to  the  incident  wave  with  the  energy  density  E, 
which  we  may  assume  to  be  perfectly  absorbed  by  the  surface.  This 


leads 


to  a  term  P  „  =  E.  The  other  part  is  due  to  a  reflected 
rad  1 


2  =: 


as  re- 


wave  with  the  energy  density  ^  E,  which  we  may  imagine 

—  p  — 

emitted  by  the  surface;  this  causes  a  pressure  Pracj  2  ~  ^  ji* 

The  whole  pressure  P  =  (t  +  V^)  E  may  therefore  be  regarded 

Feld.  Ij  * " 

as  if  the  surface  struck  by  the  plane  wave  absorbed  its  energy  per- 
fectly  and  reenitted  the  amount  ^  *  E  (see  ref.  5,  p.  294) • 

This  statement  gives  of  course  not  a  description  of  what  actually 
happens.  A  physical  explanation  of  the  actual  process  at  the 


reflector  is  offered  in  Sec.  12  below. 

c.  Radiation  pressure  at  oblique  incidence. 

If  the  plane  wave  strikes  a  wall  inclined  at  an  angle 
tS  with  respect  to  the  direction  x  of  the  wave  propagation  (Fig.  49 
Page  7)  the  radiation  pressure  upon  this  wall  is  found  by  trans¬ 
forming  the  fixed  coordinate  system  x  y  z  into  the  system  x’  y’  z 
(Fig.  5,  Page  8)  whose  x’y’  axes  are  rotated  through  the  angle  ziL 
Applying  the  transformation  formula  (1.13)  for  the  tensor  T,  we 
get  the  new  tensor  components  in  the  x1  y!  z  system.  From  (1.1 3) 
and  (2.49)  we  find  for  a  finite  bean  traversing  a  liquid  with  ' 
static  pressure  pQ  and  striking  a  -perfectly  absorbing  surfaces 


rp  |  _ 


-(p0  +  E  cos^"$) 

|  sin  2  ijl 

O 

0 

1  sin  2  tX 

2 

-(pQ  +  E  sin~  £-$) 

0 

0 

0 

V, 

(2. 


Fig.  10  shows  the  stresses  exerted  upon  a  volume  element 


oriented  according  to  the  new  axes  x*  y*  z0 


y'  v 


pn  +  E  sin^-S*  \  \ 


\  ^O 

\  •« 


pQ  +  E  cos 

\ 

\  !sSs  2A\  /  v\  ^ 

'  \v  >4  \  \ 

>V"  \  v 

N  \ 

p  +  E  cos^X  \ 

0  \  \ 


2,-a 


>x 


pD  +  E  sin^gP 


Fig.  10 

Components  of  stress  due  to  radiation  at 
a  volume  element  inclined  at  an  angle  iSl 
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If  we  imagine  a  naterial  surface  inserted  in  the  way  of  the  beam 
at  an  angle  1$  to  be  adjacent  to  the  volume  element  in  Fig,  10  and 
parallel  to  the  y'-axis,  it  undergoes  a  pressure  pQ  +  E  cos2  2$  in 
the  direction  of  its  normal,  x',  since  it  has  to  exert  the  same 
pressure  upon  the  y1  z  surface  of  the  volume  element.  The  radia¬ 
tion  pressure  normal  to  the  perfectly  absorbing  surface  is  there¬ 
fore  given  by  E  cos2  2%  If  we  assume  a  reflecting  surface,  we 
have  to  add  the  radiation  pressure  exerted  by  the  reflected  wave 
regarded  as  emitted  in  the  direction  shown  in  Fig.  11, 


Reflected  have 


Fig.  Ill 

Direction  of  propagation  of  the  reflected  wave 
for  oblique  incidence 


As  can  easily  be  seen,  we  have  only  to  change  In  Eq.  (2.56)  into 
-  2$  to  get  the  stress  components  related  to  the  reflected  ■'.rave. 

The  pressure  exerted  by  the  reflected  wave  normal  to  the  reflecting 


surface  is  therefore  ^  **  E  cos""  .  The  total  radiation  pressure 


/  2  ” _ 2 

normal  to  the  reflecting  surface  amounts  to 


(U  =(1+-j'2)Ecos22fi 

x1 


(2.57) 
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As  Sq.  (2.56)  shows,  the  wave  exerts  also  a  shearing  force  upon 
the  reflector,  if  the  latter  is  inclined  at  an  angle  2$  ?  it 
amounts  to  ~  sin  2  4$  and  therefore  has  its  maximum  at  2$.  =  45°. 

Owing  to  the  reversal  in  the  sign  of  for  the 
reflected  wave,  the  shearing  force  due  to  this  wave  is 
-  Y  2  |  sin  2  2$  ,  The  resultant  shearing  force  per  square  centi¬ 
meter  connected  with  the  normal  radiation  pressure  in  (2.57)  is 
therefore 

s  _ 

O'  -  y2)  ?  Sin  2  zfi  (2.57a) 


,'s  \  = 

v'  ra<Vy, 


10 


It  vanishes  for  a  perfect  reflector  (  y  =  1 )  or  for  2$.  =  0. 

It  may  be  noted  that  the  reflection  coefficient  \/  is 
a  complicated  function  of  the  angle  of  incidence  2$-,  because  a 
compressional  wave  striking  the  wall  at  the  angle  induces  in 
general  compressional  and  rotational  waves  inside  the  reflector. 

Our  object  here  is  limited  to  the  knowledge  of  the  average  value  in 
time  of  the  resultant  second  order  forces  acting  upon  the  reflector 
and  due  to  radiation  pressure;  the  reflector  is  assumed  to  be 
characterized  by  the  amplitude  reflection  coefficient  y  (  2&),  and 
to  absorb  all  the  energy  not  reflected,  so  that  no  radiation  leaves 
the  rear  of  the  reflector.  In  the  special  case  of  a  perfect  re¬ 
flector  \/  —  Ti  and  (by  definition)  is  independent  of  . 

In  the  case  where  some  of  the  energy  is  transmitted 
through  the  reflector,  we  can  assume  that  the  reflector  radiates 
this  part  from  its  rear  face.  The  reactional  force  is  proportional 


to  the  energy  passing  through  the  reflector,  and  of  opposite  sign 
to  the  forces  of  both  the  incident  and  the  reflected  wave  at  the 
front  side.  If  xie  call  o'  the  amplitude  transmission-coefficient 


of  the  wave  leaving  the  rear  of  the  surfaces,  the  resulting  radia¬ 
tion  pressure  is  obviously  given  by 


5rad  =  (1  +  ^  2  -  6  2)  E  (2. 

for  a  wave  at  normal  incidence.  For  oblique  obstacles  the  corres¬ 
ponding  forces  can  be  found  in  the  same  way. 

As  an  example  of  the  action  of  the  radiation  pressure 
upon  oblique  surfaces,  we  compute  the  force  exerted  upon  the  device 
shown  in  Fig.  12,  which  is  used  frequently  for  measuring  radiation 
pressure. 


Normal  (li)  and  shearing  (S)  forces  at  a 
wedge  undergoing  radiation  pressure 


The  width  of  the  vane  in  Fig.  12  is  X.  ?  the  dimension  perpendicular 
to  the  plane  of  the  paper  is  JL-,  the  reflection  coefficient 
The  whole  force  normal  to  one  side  of  the  vane  is  II  =  JX  (I  +  y2)  E 
ccs^  2$,  the  shearing  force  is  S  =^CA(li  -  N^)  ~~  s^n  2  (We 

neglect  the  diffraction  effect  induced  by  the  edges  of  the  vane  and 
effects  due  to  the  hydrodynamic  flow).  As  we  are  interested  in 


54 


the  resultant  force  in  the  direction  x  of  wave  propagations  we 
find  (Fig.-  12): 


P  -  2  S  cos  (90  -%$)  t  2  II  cos  2$  =  2(S  sin  zJ?  +  N  cos  2$  ) 

=  2  X${(1  -^2)  §  sin  2  z8  sin  Z&+  (1  +^2)  I  cos2 2?  cos  aPj 

=  2: ijii  cos^(l+y2  cos  2  zS ) 

«  '  ■ 

Introducing  the  angle  a  —  2.(90°  -  zX)  —  it  -  2  2a included 
between  the  two  sides  of  the  vanes  and  calling  attention  to  the 
fact  that  v  is  a.  function  of  a  by  writing  V  (a)  instead  of  \7  s 


V 


we  have: 


_  ^  O  * 

P  =  2v(Ae  •  sin  ^  (l  -  v  (a)  *  cos  a) 


Y 


(2.59) 


a  =  180°  corresponds  to  a  plane  surface  normal  to  the  x-axis  and 
from  (2.58)  we  find  as  expected*  P^.  =  2.J/JiE  (1  +  y^2)*  The  i°rce 
at  the  vane  tends  toward  zero  with  decreasing  a. 

It  might  be  interesting  to  note  that  for  a  =  90° 


P  becomes  independent  of  y'  (tt/2.)  and  equals  (Px) 


6 


x'90° 


\ fT  Xfj. 


For  every  value  oi 


a  90°-vane  undergoes  the  same  force  (P  )  as 


Y  ’ "  '■  - ° x"x' 

a  perfect  absorber  presenting  an  area  I f~2.  X normal  to  the  incident 
wave.  Since  the  reflected  wave  propagates  in  this  case  in  the 
direction  normal  to  the  x-axis*  it  carries  no  mean  flux  of  momentum 
along  the  x-direction.  This  explains  the  independence  of  (Px)gQ0 
from  the  reflection  coefficient  ^  (7r/2), 


o,  Fnaiation  Pressure  on  a  Perfect  Absorber 
at  Finite  Amplitudes 

All  previous  considerations  are  valid  only  for  small 
amplitudes*  where  k^o<<1i,  because  we  used  a  development  in  series  of 
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the  exact  solution  for  u  and  p,  retaining  terms  up  to  the  second 
order  only.  But  since  we  have  the  exact  solutions  for  liquids 
in  Eqs.  (2.30)  and  (2.31 )s  we.  are  able  to  deal  also  with  finite 
amplitudes.  lie  limit  our  consideration  of  finite  amplitudes  here 
to  the  case  of  a  wave  traveling  in  one  direction  only's  terminating 
at  a  perfect  absorber.  This  example  will  show  clearly  enough 
the  difference  between  the  physical  characteristics  of  small  and 
finite  amplitudes. 

For  computing  the  tensor  components  in  Eq.  (2.18)  we 

•—  O 

need  the  mean  pressure  p  and  the  mean  flux  of  momentum  ^>u  . 

Using  the  solutions  (2.30)  and  (2.31)  for  u  and  p  with  y  ~  0S 
and  the  density^©  from  (2.28),  we  find 
£ 

cos  (cot  -  kx  +  lc£0  sin(cot  -  kx)  )  dt 


-.1  ( 
P  =  T 


P  Jr 


2  _  1 


fu 


■p 


‘ofi  cos 2  (cot  -•  Inc  +  k^n  sin  (cot  kx)  ) 

1  -  k^0  cos  (cot  -  lac  +  kq0  sin  (cot  -  kx)  ) 


dt 


n  2  / 

Introducing  <xr ^>Q  ~  k  /p  and  the  new  variable  cot  -  lex  =  "2r,  we  have 


“  _  1 

P  p  *  2tt 


cos  ("£  +  k^Q  sin  )  dTf 


(2.60) 


k2^o2  #  1  ^  cog2  (r  +  sin  -V  ) 

P  2ir  yQ  1  -  k£o  cos  ( X  +  k^Q  sin  'If  ) 


(2.61) 


The  integral  in  (2,60)  is  found  to  be  represented  by  the  Bessel- 
f unction  of  index  1,  -J  (k^),  So  we  have 

5 "  -  y  Ji  «.> 


(2.62) 
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For  small  values  of  (k^<Kl )  the  function  is  given  by 


j  (kF  )  =  -42  - 
1  so'  2 


Ko  C*o>3 


16 


+  «... 


Taking  the  first  term,  we  get  for  small  amplitudes  from  (2,62) 


according  to  Eq.  (2.42).  For  larger  amplitudes  the  mean  excess 
pressure  p  deviates  from  —  E,  since  p  is  represented  by  (2.62). 
We  find: 

k?Q  =  0.05  0.1  0.3  0.5  0.7  0.8  0.9  (1.0) 

-§p  =  O.OOt!  C.005'  0.044  0.12T  0.230  0.295  0.365  (0.440) 


(A) 


As  k^Q  approaches  1 ,  the  values  given  here  for  -  ppo  are  only  of 
theoretical  interest  since  the  concept  of  constant  compressibility 
loses  its  meaning  at  such  high  compressions  as  are  encountered  at 
larger  values  of  k£0« 

The  expression  (2, 6n)  for  P u^  =  2  cannot  be 

represented  by  known  functions  in  a  closed  simple  form,  as  is 
the  case  with  p. 

The  integral 


•  1  r*  cos?(?-t  ga  stof) _ 

J-1  (k£Q)  “  Sr  jQ  1  -  k^Q  cos(r+  k^Q  sin^  )  d 


(2.63) 


has  been  evaluated  numerically  for  different  values  of  the  para¬ 
meter  k£  .*  The  result  is  as  follows: 


*  The  writer  is  grateful  to  Dr.  S.  P»  Frankel  of  this  institution 
for  suggesting  a  quick  and  simple  method  for  evaluating  the 
integral  in  Eq.  (2.63). 


57 


k?0  =  0.05  0.1;  0.3  0.5  0.7  0.8  0.9  1i.O 

F(k?0)  =  0.500  0.5025  0.524  0.571  0.685  0.786  11.057  CO 


(S' 


For  small  values  of  the  integral  F(k^0)  equals  l/2  and 

p)u2  =  k2£o2/2p  =  E  according  to  Eq.  (2.61).  For  increasing  values 

of  the  value  of  pu2  undergoes  a  deviation  from  the  energy 

density  E.  lie  obtain  for  p«  ^>u2  =  (k?0)2, 

kl  =  0.05  0.1  0.3  0.5  0.7  0.8  0.9  C1*0) 

2 _  (C 

p»  pu2  =  0.001  0.005  0.047  0.143  0.336  0.503  0.856  (°0 


The  quantity  p  +  ^)u2  changes  thus  with  increasing  amplitudes: 
kl  =  0.05  0.1  0.3  0.5  0.7  0.8  0.9  (l.O) 

p*  (p+ pu2)=  0.000  0.000  0.003  0.021  0.106  0.208  0.491  (<*>) 


(d; 


For  smaller  values  of  k£0  the  following  series  developments  may 
be  used: 


pp  =  -  j|(k?o)2 
p^>u2  =  +  ^(k?o)2 


+  16"^o)4  “  384(k^o)6  +  T8432(lc^o)S”" 

+  l(k?  )4  +  %1: K  )6  +  -^|(k8  )S  +  .. 
4  o'  96  °  3072  0 


(2.62 

(2.61a 


The  development  of  pp  follows  immediately  from  (2.62)  and  the 
well-known  series  development  of  the  Bessel-function  J^(x).  The 
development  of  p  pu2  was  obtained  by  developing  the  integrant  in 
(2.61)  in  powers  of  k^Q  cos(^+  kZ,Q  sin?),  converting  the  powers 
of  the  cosines  into  cosines  of  multiples  of  the  argument,  then 
integrating  term  by  term  by  using  the  relation 

(l/2ir)  cos  n(T  +  kE  sintr)  dV  —  (-ti)  Jn(nk?0) , 

°o 

and  finally  developing  the  Bess el-functions  Jn(nk£o)  in  powers 
of  nk£Q.  The  physical  conclusions  to  be  drawn  from  these  re¬ 


sults  are  as  follows: 
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a)  Infinitely  extended  plane  comoressional  waves, 
not  communi eating  with  undisturbed  regions 

For  small  amplitudes,  as  table  (D)  shows  us  again,  the 
two  terms  p  and  pu2  cancel  each  other,  because  both  quantities 
equal  the  mean  energy-density  -  E  or  +  E.  At  increasing  ampli¬ 
tudes  the  sun  (p  +  pu2),  which  according  to  (2,1 8)  is  responsible 
for  the  radiation  pressure  in  this  case,  is  positive.  At  higher 
amplitudes,  therefore,  even  in  the  case  of  an  infinitely  extended 
plane  wave  in  a  non-absorbing  medium,  a  radiation  pressure  exists, 
though  its  value  is  negligible  if  the  amplitudes  are  not  extremely 
high.  In  the  limit,  at  which  k^Q-y  1 ,  this  radiation  pressure 
would  theoretically  tend  towards  infinity.  This  can  be  understood 
physically  by  the  fact  that  for  =  1  the  density  becomes 
infinite  once  in  each  cycle,  leading  to  infinite  values  of  the 
double  mean  kinetic  energy-density  pu2.  But  at  such  great  ampli¬ 
tudes  the  concept  of  a  constant  compressibility  will  surely  not 
hold.  Therefore  these  statements  concerning  extremely  large 
amplitudes  (k£0-»l)  are  only  of  theoretical  interest. 

b)  Finite  plane  beam  surrounded  by  or  communicating 
with  undisturbed  regions  and  incident  on  a  perfect 
absorber 

O 

In  this  case,  as  stated  above,  only  the  term  p u^  is 
responsible  for  the  radiation  pressure,  as  the  change  in  pressure 
p  will  be  equalized  by  the  hydrostatic  pressure  pQ  outside  the 
finite  beam.  For  small  amplitudes,  as  may  be  seen  from  Eq.  (2.6l) 
and  the  foregoing  table  (C),  and  as  is  known  from  previous  considera¬ 
tions,  pu2  =  E.  That  is,  as  a  first  approximation,  the  radiation 
pressure  equals  the  mean  total  energy  density  E.  At  higher  ampli- 


tudes  this  is  no  longer  true. 


The  general  statement,  correct  for  any  amplitude,  in 
this  case  is  as  follows: 

The  radiation  pressure  exerted  by  a  compressional  plane 
wave  which  is  surrounded  by  or  in  communication  with*  a  medium 
not  affected  by  the  wave  motion  and  in  contact  with  the  rear  of 
a  totally  absorbing  surface ,  equals  2  when  is  the  mean 
kinetic  energy-density  of  the  plane  wave  in  Eulerian  coordinates. 
This  is  correct  for  any  amplitude,  since  u^  =  2  E^-j^  by 
definition.  With  this  in  mind,  the  formula  already  given  in 
Eq.  (1.17)  can  be  regarded  as  generally  correct. 

As  to  mean  energy  densities  at  finite  amplitudes,  one 
must  bear  in  mind  that  E  ^  and  E^.^  and  also  E|0^  and  are 

not  the  same  in  Eulerian  as  in  Lagrangian  coordinates.  This 
disparity  disappears  at  small  amplitudes,  that  is,  in  the  second 
order  of  approximation.  The  acoustic  intensity  is  given  by 
J  =  p^u^,  using  our  previous  consideration  concerning  a  perfect 
absorber  (page  4l)«  The  mean  intensity  at  any  amplitude  is 
therefore  given  by 


*  Such  a  communication  can  be  regarded  as  realized  for  instance 
by  a  small  hole  in  the  absorbing  surface,  on  the  front  of  which 
is  the  incident  beam,  while  the  rear  face  is  in  contact  with 
an  undisturbed  region  of  static  pressure  p  .  As  G.  Richter"1 ^ 
has  already  correctly  remarked,  the  hole  sfiould  be  visualized 
as  opaque  to  acoustic  radiation  but  as  allowing  the  equalization 
of  average  pressure  between  the  two  regions,  as  might  be  effected 
by  having  a  small  totally  absorbing  piston  free  to  move  in  the 
hole. 
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-  T  r  2  2 r 

J  =  1  V  dt  =  p2-  *  ~  \  cos2  (cut  -  kx)  dcot 

"p  vo  Jo 

.2'  2  - 

—  k  o  ~~  ? 

or  J  =  c  •  =  c  o  u*^  =  2  c  Es 

2£  \  0  pot 

•using  Eqs*  (2.24)?  (2.25)?  (2.40)  and  c ^§pQ  -  1.  The  amplitudes 
belonging  to  a  mean  intensity  J  are  therefore  always  given  by 


icr  =S/3_|_S 


If  we  as  suit. e  radiation  in  water  at  20°  C?  with  k£  =1 

o 

(a  purely  hypothetical  value)?  the  intensity  J  is  found  to  be 


1.49(1 0‘5)  cm  kg 
2  v 46(10“”)  sec. cm2 


1.62(109) 


9.81  watt 
1 00  cm^ 


or  J  =  1.59(10^) 

This  expression  shows  how  far  we  always  are  experimentally  from 

-2 

"finite"  amplitudes,  because  even  an  intensity  of  100  watt  cm 
is  to  be  regarded  as  very  high  in  water.  The  value  of  klQ  for 

an  intensity  J  =100  watt  cm“^  is  found  to  be  7.93(10-^). 

For  all  practical  purposes  in  liquids  the  amplitude  as 
well  as  the  quantity  klQ  can  be  regarded  as  small.  The  radiation 
pressure  is  therefore  described  for  all  cases  which  are  experi¬ 
mentally  feasible  by  the  considerations  in  Section  7?  where  we 
dealt  with  small  amplitudes. 

Nevertheless  the  theoretical  results  for  finite  ampli¬ 
tudes  are  not  without  interest?  since  they  show  that  the  relation 
P_^d  =  mean  total  energy  density  E  cannot  be  regarded  as  a  "basic" 
physical  law  as  it  is  considered  to  be  in  electrodynamics.  For 
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mechanical  wave  motions  a  law  of  this  kind  holds  only  at  small 
amplitudes,  whereas  in  general  at  all  amplitudes,  the  radiation 
pressure  is  linked  with  the  expression  for  the  mean  density  of 
kinetic  energy  in  Euler ian  coordinates. 


As  E  =j2k.n  4  Epot,  the  difference  E  -  2%^  =  Epot  - 
E^.^  =  1/2.  (  <oq  -  p u^)  in  agreement  with  Sa.  (2.40).  In 
liquids  this  difference  turns  out  to  be  small  of  the  sixth  order 
in  (kqQ).  In  order  to  show  this,  one  uses  the  series  development 
of  p  pu2  as  given  by  Eq.  (2,62)|  £>  p  u^  is  found  with  the  aid  of 
(2.30)  to  be  represented  by  l/2  (k21)2  1:  4  J  (2k£  )  j.  Making 

u  O  j 

use  of  tlie  known  series  development  of  Jo(x)  =  x^/2  “  ky24  +  •  •••, 
we .  obtain  indeed  j5  (E  -  22^..^)  =  -  (k^D)^/64  4  ....,  that  is, 
only  terras  of  the  sixth  order  and  higher  are  involved. 


9.  Note  on  an  Integral  of  the  Eulerian  Equation  of  Motion 
and  the  Flow  of  Mass  for  Infinitely  Extended  Plane 
Haves  in  Liquids 

For  liquids  with  constant  compressibility  we  found  that 
the  solutions  of  the  Lagrangian  equation  (2.22)  and  consequently 
those  of  the  Eulerian  equation  (2.14)  were  functions  of  the  argu¬ 
ment  (cot  4  kx)  only.  Using  the  fact  that  therefore 
3/dt  =  4  C  •  c>/dx,  Eq.  (2.14)  can  be  integrated  with  respect  to 
x  or  t  and  vie  obtain 

4c*  ^>u  4  ^)U^  4  p  =  C  (2.64) 

C  is  independent  of  x  and  t,  but  not  necessarily  of  the  amplitude 
of  the  wave  motion.  He  can  prove,  however,  that  C  has  the  value 
zero.  This  proof,  which  we  will  now  give,  verifies  at  the  same 


time  the  fact,  already  used  in  Part  II,  Sec.  5,  that  any  function 
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f(cct  T  x),  which  was  found  to  satisfy  the  Lagrangian  equation 
(2.22),  leads  to  an  exact  solution  of  the  corresponding  Eulerian 
equation  (2.14) $  when  the  transformations  of  the  quantities 
p*5  P*  into  the  corresponding  Eulerian  quantities  us  p,  p 
are  properly  made. 

Let  32)/  ^  t  =  u*  =  u*  (cot  +  lac)  be  the  general  solution 
in  the  Lagrangian  system.*  Prom  Eqs.  (2,10)  and  (2.28)  we  obtain 


P* 


13?  1  3£ 

“  p  sx  ~  ±  pc  at  -  ±  f0c 


u* 


(2.65a) 


P  o 


(2.65b) 


Transforming  p#,  and  u*  into  p,  and  u  by  replacing  the  variable 
x  in  u*  by  (x  —  O  and  inserting  p  and  p,  as  found  in  Eqs.  (2.65&sb)_, 
into  Eq.  (2.64) 5  we  have 


T  c 


-E° 


U 

4_a_ 


+  H 


1  + 


_  U  i  foc 


u  - 


c 


or  C  =  0 


Thus s  Eq.  (2.64)  provides  us  with  the  general  relations  valid  in 
liquids  at  any  tine  t  and  at  any  x: 


p  + 


+  c  *c>u 


(2.66) 


Averaging  (2.66)  in  times  the  mean  mass  flow  carried  by  the  wave 


u*(cot  +  kx)  stands  for  u*(oot  -  kx)  +  Up  (cut  +  kx)s  where  u.  ^ 
are  arbitrary  functions  of  the  arguments  (cot  +  lac).  5 
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motion  through  a  fixed  unit  cross-section  turns  out  to  be 

o  u  =  +  ^  (p  +  pu^)  (2.66a) 

The  plus  or  minus  sign  corresponds  to  waves  traveling  to  the  right 
or  left.  The  value  of  (p  +  pu2)  at  finite  amplitudes  has  been 
given  in  Table  (D),  p.  57.  From  Eqs.  (2.60)  and  (2.6l)  we  obtain 
in  the  case  of  sinusoidal  wave  motions 


1 


r 


>2tt 


-  (p  +  pu2)  =  . 

c  y  c 


\ 


COS  (F+  +  k?n  sin  77) 


1  +  k l  cos  (f  +  k?  sin  77  ) 

u  •—  O 

The  series  development  of  the  integral  in  powers  of  (k^o),  is 
found  by  adding  Eqs.  (2.62a)  and  (2.61a)  on  page  57: 

■p  =  ±  ~  (P  +  £>u2)=  +  f>Qc  (^  (k?0)4  +  }££  d^0)6  +  ...) 


dt  (2.67) 


(2.68) 


Eq.  (2.68)  brings  once  more  to  evidence  the  result  formerly  ob¬ 
tained,  that  for  small  amplitudes,  that  is  disregarding  terms  of 
higher  than  second  power,  the  sum  p  +  pu2  cancels  out.  It 
confirms  also  the  conclusion  made  concerning  the  "Constant"  in 
Section  7b,  namely  that  this  constant  is  independent  of  (ksQ)  to 
the  second  order  of  approximation. 

The  series  development  of  the  change  in  average  density 
due  to  the  wave  notion  p=  p  “  p0  follows  likewise  immediately 
from  (2.68).  Eq.  (2,65b)  shows  that  p*  -  -  +  p*u*/ c . 

Transforming  this  into  Eulerian  coordinates  and  talcing  the  time 
average  we  obtain 


Af  =  fPo=±  f u/o 


fo  +  Ss  (k?0>6  +  •••) 


(2.68a) 
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As  stated  on  page  44  the  change  in  mean  density  41  p  vanishes 
to  the  second-order  approximation;  it  is  of  the  fourth  order 
in  (ksQ)  and.  positive,  whereas  the  corresponding  change  in  pres¬ 
sure  p  -  -  E  is  of  the  second  order  ana 

The  average  flow  of  mass  p  u  vanishes  to  the  second 
order  of  approximation,  as  stated  in  various  papers  (see  for 
instance  reference  16),  This  does  not  mean,  however,  that  it  is 
strictly  zero  as  is  seen  from  Eq.  (2.68)*.  In  liquids,  the 
values  of  k£  that  can  be  attained  experimentally  are  very  small, 
and  the  quantity  pu  amounts  numerically  to  an  extremely  small 
value.  It  may  be  noted  that  the  existence  of  a  "net  flow  of  mass11 
does  not  necessarily  imply  an  actual  transport  of  matter ,  that  is, 
a  mean  particle  displacement.  Indeed,  in  a  pure  sinusoidal 
wave-motion  every  particle  is  found  at  exactly  the  same  place  as 
before  after  a  full  period.  The  average  values  in  time  of  £  or 
u*  are  zero  in  this  case.  Nevertheless,  owing  to  the  fact  that 
p  is  not  a  constant,  but  ~  P2  +  **‘,!’  ^1s  Po**’ 

denoting  higher-order  terms  in  p  ,  varying  sinusoidally  in  time) 
the  time-average  of  the  "flux  of  mass -density”,  given  by  (^  *u), 
may  have  a  value  different  from  zero. 

The  quantity  p  u  is  independent  of  x.  To  prove  this, 
we  take  the  time-average  of  the  Eulerian  equation  (2.14)  and  find 
for  a  wave  motion  periodic  in  time, 


(p+  (p+  r2)  =  0 


*  The  only  exception,  where  cu  equals  zero  strictly,  is  the  case 
of  a  standing  wave.  ' 
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or,  according  to  Eq.  (2.66a),  ^ 

pu  is  independent  of  x  and  a  function 

(l£o)  only- 


0.  Consequently 
the  wave  cutipi-itude 


thus  we  have  in  liquids  at  any  amplitude  for  an  in- 
£Ar,4+<ily  extended  plane  wave 


with 


c  • 


f(lc?D) 


(2.69) 


w4^  +  5^,0  +  i2^ . 


17 

728 


4567 


,8 


10.  Radiation  Pressure  and  Viscous  Absorption 
The  Eulerian  differential  equation  (2.14)  includes 
the  case  of  plane  compress ional  waves  in  viscous  fluids.  But  in 
the  case  of  viscous  absorption  an  additional  term,  associated  with 
the  coefficient  of  viscosity  no,  enters  the  differential  equation 
both  in  the  Eulerian  form  Eq.  (2.14)  and  the  Lagrangian  form  Eq. 
(2.22).  D’Alembert's  general  solution  f(cot  T  kx)  is  then  no 
longer  a  solution  of  Eq.  (2.22),  as  is  well  known,  and  we  cannot 
make  use  of  the  relation  c)/  <)t  =  +  c  -~6/c)x.  The  conclusions 
based  on  this  relation  in  the  foregoing  section  9  therefore  do  not 
hold  for  viscous  liquids.  Still,  talcing  the  average  value  in  time 
of  Eq.  (2.14),  as  we  did  at  the  end  of  section  9,  we  find  also  in 
the  case  of  viscous  liquids  (or  gases)  that  B  (p  +  pu2)/<3x  =  0 
for  a  wave  notion  periodic  in  time  and  therefore  p  +  pu^  =  constant. 
As  the  amplitude  of  a  uni&irect ional  plane  wave  decreases  steadily 
along  its  direction  of  propagation  through  a  viscous  medium,  the 
values  of  u  and  p  steadily  approach  zero.  Thus  the  value  of  the 
constant  is  necessarily  zero  and  p  and  pu2  cancel  each  other  in 
the  same  way  as  they  do  in  absence  of  viscosity.  The  radiation 
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pressure  p  +  p  due  to  a  beam  of  infinite  width  (and  of  any 
amplitude)  is  therefore  zero  at  a  perfect  absorber  in  the  case 
of  a  viscous  liquid  (or  fluid).  The  same  conclusion  holds  for 
a  reflected  plane  i/ave  leaving  an  absorber. 

For  a  beam  of  finite  width  the  problem  becomes  more 
involved.,  owing  to  the  interaction  of  the  beam  with  undisturbed 
regions.  For  non-viscous  liquids ,  the  time-average  pressure  p 
is  constant  along  the  beam,  or,  even  when  periodic  in  space, 
still  its  average  value  along  the  beam  is  constant.  This  fact 
has  enabled  us  to  deal  with  this  interaction  in  a  relatively 
simple  manner. 

On  the  other  hand  in  viscous  liquids,  the  value  of  p 
in  the  beam  changes  exponentially  along  the  beam,  and  the  inter¬ 
action  with  the  undisturbed  region  outside  the  beam  involves 
rotational  motion  of  the  fluid.  The  problem  of  radiation  pres¬ 
sure  becomes  closely  linked  with  that  of  the  hydrodynamic  flow. 

A  very  complete  treatment  of  the  forces  due  to  second  order 
effects  in  viscous  media  is  given  in  the  fundamental  paper  by 
Eckart. ° 

11,  Rayleigh  Pressure  and  Langevin  Pressure  in 
Liquids 

In  papers  dealing  with  acoustic  radiation  pressure 
special  expressions  have  recently  been  introduced  in  the  litera¬ 
ture  to  denote  the  special  circumstances  under  which  radiation 
pressure  may  be  observed.  It  is  our  purpose  to  clarify  the 
physical  meaning  of  these  usages.  According  to  our  previous 


considerations  we  can  distinguish  four  different  cases  of 
pressures:  they  are  connected  on  the  one  hand  with  the  coordinate 
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system  used  (Eulerian  or  Lagrangian) »  and  on  the  other  hand  with 
the  kind  of  interaction  of  the  acoustic  beam  with  an  undisturbed 
mediun9  that  iss  whether  the  beam  is  regarded  as  of  infinite 
width  and  not  communicating  with  an  undisturbed  part  of  the  mediums 
or  of  finite  cross  section  and  surrounded  by  the  undisturbed 
medium.  lie  may  clarify  these  cases  as  follows: 

1.  (L.  -  i. ):  Lagrangian  coordinates  and  infinite  plane  wave. 

2.  (L.  -  c.):  Lagrangian  coordinates  and  wave  region  communicating 

with  undisturbed  medium. 

3.  (E.  -  i.):  Eulerian  coordinates  and  infinite  plane  wave. 

4.  (E.  -  c.):  Eulerian  coordinates  and  wave  region  communica.ting 

with  undisturbed  medium. 

We  will  deal  with  these  cases  successively  (always  assuming  normal 
incidence): 

1.  (L.  -  i. )  At  any  point  on  a  surface  following  the 
motion  of  a  particle  in  a  liquid 9  the  pressure  p#  varies  purely 
sinusoidally  (Sq.  (2.25)  )•  Its  time  average  value  is  therefore 
zero:  p#  =  0.  An  observer  moving  together  with  a  particle  would 
register  a  mean  pressure  equal  to  zero.  The  mean  pressure 

in  the  medium  for  a  stationary  observer  would  be  found  to  be 

lowered  to  pQ  —  E  (1  +  y^)  at  small  amplitudes  in  an  infinitely 

extended  bean,  where  p  is  the  pressure  in  absence  of  wave  motion, 

o 

For  larger  amplitudes  the  decrease  in  mean  pressure  follows  from 

Eq.  (2.62). 

2.  (L,  -  c. )  The  mean  pressure  p*  for  a  moving  observer 
is  zeros  just  as  in  1 ,  (L.  -  i. ).  For  a  stationary  observer  the 
mean  pressure  p^=  pQ,  since  the  lowering  of  pQ  by  the  wave  notion 
is  now  counter-balanced  by  the  action  of  the  medium  of  static 
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pressure  pQ  which  surrounds  the  beam  region. 

3.  (E.  -  i.)  The  mean  pressure  jL  at  fixed  coordinates 
is  lowered,  as  already  mentioned  under  1.  (L.  -  i. ) ,  by  the 
amount  —  E  (l  +  y2)  at  small  amplitudes.  If  its  value  was  p0 
in  absence  of  wave  motion,  then  for  an  infinitely  extended  plane 
wave,  not  communicating  with  any  undisturbed  medium,  the  value  is 
pQ  ~-E  (1  +  y2).  The  radiation  pressure  Pra(j,  which  we  have  been 
led  to  identify  with  (p  +  ^>u^)  in  Eulerian  coordinates,  becomes 
zero  for  small  amplitudes.  For  larger  amplitudes  numerical  values 
of  the  radiation  pressure  can  be  found  from  Table  (D)  on  page  57 
or  from  Eq.  (2.69). 

4«  (S.  -  c. )  If  the  acoustic  beam  communicates  with  a 
medium  not  affected  by  the  wave  motion,  the  mean  pressure  p-j. 
becomes  equalized  to  the  value  of  the  static  pressure  pQ  of  the 
undisturbed  medium.  This  raises  the  mean  density  in  the  beam 
by  -A.jp’  =  E/cj2  (l  +  Y2)  at  small  amplitudes,  as  follows  from 
Eq.  (2.48b).  The  radiation  pressure  =  E  (1  +  y2)  for 

small  amplitudes. 

The  notation  "Rayleigh  pressure"''  (Ppayqeigh^ s  as  use<^ 
in  the  literature  8*1} 15,  means  the  average  excess-pressure,  due 
to  the  wave  motion,  which  would  be  noticed  by  an  observer  moving 
with  a  particle.  It  is  therefore  identical  with  the  quantity 
p*  in  our  notations 


P  —  p* 

Rayleigh  t 


(2.70) 


For  liquids  with  constant  compressibility,  p*  and  consequently 
the  "Rayleigh-pressure"  is  zero  at  all  amplitudes.  For  media 
that  have  a  more  complicated  relation  P*(p)  between  pressure 
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and  density  -  as  in  gases  for  instance  -  p*  and  therefore  the 
Rayleigh-pressure  is  different  from  zero.  (See  Part  III,  Sec.  3  and  5) 
The  expression  "Langevin  pressure”  (PLangevin)  use(^s 
following  Hertz  and  Hende  ,  for  the  difference  between  the  pressure 
p*  observed  at  a  moving  particle  or  plane  and  the  mean  pressure  p 
at  fixed  coordinates: 


langevin  ^  ^  ^Rayleigh  “  ^ 


(2.71) 


For  a  plane  conpressional  wave  traveling  in  one  direction  in  a 
liquid,  we  found  at  small  amplitudes  p*  =  0  and  p  =  —  E,  therefore 

=  =  E  [^o  «!]  <2'72) 

This  result  turns  out  to  be  independent  of  the  nature  of  the 
medium,  that  is,  independent  of  the  special  function  p( ^> )  connec¬ 
ting  pressure  and  density.  This  fact  has  already  been  stated  in 
3,1 


various  papers 


An  exact  and  simple  proof  is  as  follows: 

.  =  _ 

“b t2 


For  small  amplitudes  we  have  according  to  Eq.  (2.21) 


p(x)  =  p*(x  -  K)  =  p *-(x)  -  K 


- 


t  0  0  0  9 


Therefore 9  at  snail  amplitudes 


P*(x)  -  p(x)  =  (?  ~ ) 


s2* 


Now  we  can  write 


—  (£  &£) 

at  ^  at' 


(^)2 

lat; 


(2.73) 
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and  taking  the  average  time-value  and  regarding  ^  as  a  function 
periodic  in  time,, 


Thus  we  have  from  Eq.  (2.73) 


p*(x)  -  p(x)  =+  f0  (ff)2  =  pou*2  =  <2-?4) 

For  a  unidirectional  wave  motion,  pu**2  equals  S  and  we  have 
^Langevin  =  P*  “  P  =  E.  On  the  other  hand,  Eq.  (2.74)  is  per¬ 
fectly-  general,  holding  whether  or  not  a  reflected  wave  is  present; 
but  in  this  case  pu*2  varies  from  point  to  point  and  we  must  take 
the  average  value  in  space  of  Eq.  (2»74)  if  we  wish  to  introduce 
the  mean  total  energy  density  E.  Doing  so,  we  get  from  (2.74) 


PT  .  =  p* 

Langevm 


(2.75) 


According  to  the  proof  leading  to  Eq.  (2.75),  this  general  result 
is  independent  indeed  of  the  special  function  p(p). 

In  the  literature  dealing  with  radiation  pressure  it 
is  often  the  Langevin  pressure  that  is  identified  with  the  radia¬ 
tion  pressure  exerted  upon  a  plane  obstacle  (see  references  8  and 
1).  From  the  general  result  Eq.  (2.75)  it  is  concluded  that  the 
radiation  pressure  at  snail  amplitudes  equals  the  energy  density 
in  any  fluid  medium. 

$ 

Hertz  and  Mende  seem  to  have  been  the  first  to  introduce 
the  expression  uLangevin-pressureH,  as  defined  in  (2.75).  Never¬ 
theless,  Langevin  himself  does  not  seem  to  have  had  this  quantity 
p*  -  p  (or  better  p*  —  p)  in  mind  in  his  proof  concerning  the  con- 


nection  between  iradiation  pressure  and  energy  density  as  reported 

2 

in  a  paper  by  P.  Biquard  .  According  to  Langevin's  derivations 
as  given  there s  one  is  led  to  the  conclusion  that  he  identifies 
radiation  pressure  with  p  *  -  p  *,  that  is  the  difference  between 
the  total  pressures  aspooiatej  with  noving  particles ,  at  two 
A*  <\c>-rent  regions  of  the  medium ,  one  affected  by  the  acoustic 
wave  motion  and  the  other  unaffected.  If  in  the  latter  region 
the  fluid  is  at  rest  (as  it  nay  be ,  for  instances  behind  an 
opaque  reflector),  p  **  becomes  identical  with  the  static  pressure 
Po  (see  also  the  recent  note  by  P.  J.  llestervelt17) .  In  our 
opinion  sone  conclusions  in  Biquard 1 s  report  of  Langevin's  deriva¬ 


tion  are  open  to  criticism;  moreover  we  cannot  accept  the  quantity 
pt*  ~  PQ  as  representing  the  true  radiation  pressure.  Whichsoever 
of  these  quantities  may  be  called  "Langevin-pressure",  none  of 
them  agrees  with  the  actual  radiation  pressure  given  by  (p  + 

To  diminish  the  confusion  already  widely  spreads  we  recommend 
that  the  term  "Langevin  radiation  pressure"  be  discarded  altogether. 

Acoustic  radiation  pressure  should  properly  be  identified 
with  the  expression  p  +  pu~  in  Eulerian  coordinates.  This  quan¬ 
tity  is  the  one  actually  measured;  it  is  in  fact  the  resultant 
force  per  unit  area,  due  to  the  wave  motion.  The  two  components 
p  and  jp  u  act  together;  the  contribution  of  each  depends  on  the 
characteristics  of  the  reflector  (that  is,  on  y).  The  final  sum, 
at  small  amplitudes,  is  always  E  (1  +  y2).  (See  Sec.  12) 

It  happens  that  the  value  of  Pparo.evins  as  given  by 
Eq.  (2.75),  which  is  associated  with  an  infinitely  extended  plane 


wave,  equals  numerically  the  actual  radiation  pressure  for  a  beam 
of  finite  width;  its  identification  with  the  radiation 


pressure 
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however  is  evidently  misleading  as  to  the  physical  concept  of 
this  quantity. 

12*  Lhat  xs  Measured  as  ‘‘Radiation  Pressure11  Experi¬ 
mentally? 

The  radiation  pressure  that  we  measure  experimentally 
is  physically  not  identical  with  a  ’’pressure11  in  the  hydrodynamic 
sense.  Though  of  the  same  dimensions  it  is  a  different  physical 
quantity  and  is  expressed  by  pt  -  pQ  +  p  u2  for  a  plane  compres- 
sional  wave  at  normal  incidence.  The  quantity  pu2  can  be 
interpreted  as  the  average  flux  of  mechanical  momentum  pu 
through  a  unit  area  fixed  in  space. 

Let  us  consider  a  device  D  for  measuring  the  radiation 
pressure.  It  may  contain  a  reflecting  plane  surface  struck  by 
an  incident  acoustic  beam  of  finite  cross-section.  This  ’’re¬ 
ceiving"  plane  is  to  be  regarded  as  moving  together  with  the 
immediately  adjacent  particles  of  the  medium.  This  movement, 
the  amplitude  of  which  depends  on  the  reflector,  is  purely  sinusoidal 
and  the  average  position  of  the  receiving  surface  is  identical 
with  its  position  at  rest. 

As  this  receiving  surface  is  a  part  of  the  whole  de¬ 
vice  D,  it  is  connected  mechanically  In  some  way  with  D.  The 
whole  device  D  may  be  regarded  as  having  a  large  mechanical 
Inertia  compared  with  the  inertia  of  the  receiving  surface;  its 
center  of  mass  can  therefore  be  treated  as  practically  immovable 
in  space,  even  when  the  receiving  surface  is  moving  periodically. 

For  simplicity  we  assume  that  no  acoustic  energy  leaves  from 
its  rear  side  ;  all  the  non-reflected  part  of  the  incident 
energy  is  assumed  to  be  absorbed  within  D» 


73 


a.  Ue  consider  first  the  case  of  a  perfectly  absorbing 
device  D.  The  acoustic  bean  of  finite  cross-section  is  imbedded 
in  the  undisturbed  surrounding  medium  which  also  encloses  the 
device  D,  Therefore  the  mean  pressure  is  the  sane  throughout 
the  medium,  including  the  whole  surface  of  Do  The  same  pressure 
Po  =  p+  acts  both  at  the  front  and  rear  of  D,  so  that  no  resulting 
force  due  to  a  hydrostatic  pressure  is  exerted  upon  D  in  this  case. 

Now  let  us  visualize  an  imaginary  surface  S,  enclos¬ 
ing  exactly  the  device  D  but  stationary  in  space  (Fig.  13  )• 


Fig.  13 


Schematic  diagram  of  a  device  D  for  measuring 
acoustic  radiation  pressure*  C  =  center  of  mass, 
S  =  surface  fixed  in  space  enclosing  the  whole 
device  D;  a  =  receiving  surface  struck  by  the 
acoustic  beam  and  mechanically  connected  with  D* 
p0  =  hydrostatic  pressure 


All  over  the  receiving  part  a  of  D  particles  are  crossing  the 
corresponding  part  of  S  according  to  their  periodic  movement  in¬ 
side  the  acoustic  beam.  The  whole  force  exerted  upon  the  device 
D,  and  therefore  acting  at  the  center  of  mass  C,  equals  the  gain 
in  momentum  of  D  per  second,  or  the  "flux  of  momentum1'  through  a 
per  second.  This  latter  quantity,  as  we  found,  is  given  by 
pu^  for  unit  area  (a  =  1  )•  In  the  small  time  dt  a  mass  p  u«dt 


crosses  the  unit  section  of  a  with  the  velocity  the  whole 
transport  of  momentum  through  the  section  is  consequently 
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f 


'  9 

u°dt  *  u  -  pu  -,dt« 


f 


In  unit  time  the  flux  of  momentum  amounts, 
consequently s  to  ^u^  It  may  be  noticed  that  both  the 
particles  entering  D  at  as  and  also  the  particles  leaving  D  at 
a  somewhat  later  time  contribute  to  the  force  exerted  upon  D 
in  the  same  sense  of  directions  the  departing  particles  exerting 
a  reactional  force  on  D,  According  to  our  assumption  that  no 
radiation  leaves  the  rear  of  Ds  no  force  is  exerted  at  this  side 
upon  D.  The  whole  mean  force  per  unit  area  of  Ds  struck  by  the 


compressional  wave,  therefore  equals  the  time-average  value  jO  u  , 
calculated  for  a  cross-section  fixed  in  space.  This  leads  us 
back  to  the  result  of  our  former  considerations,  that  the  radiation 
pressure  at  a  perfect  absorber  is  given  by  the  expression  ^  u^ 
in  Eulerian  coordinates.  This  double  mean  kinetic  energy  density 


f 


equals  at  small  amplitudes  the  mean  total  energy  density  E$ 
it  deviates  from  E  with  increasing  amplitude.  The  amplitudes 
that  can  be  experimentally  generated  in  liquids  are  always  "small” 
amplitudes*,  therefore  the  statement:  radiation  pressure  =  mean 
energy  density  holds  airways  for  this  case  in  liquids.  For  other 
media,  as  gases,  it  is  also  correct  for  small  amplitudes,  as  sliox-m 
in  Part  HI,  Sec.  2. 

b.  liow  let  us  consider  the  case  in  which  the  device  D 
is  not  a  perfect  absorber,  but  reflects  partialljr  or  totally  the 
incident  wave  and  absorbs  all  that  Is  not  reflected.  Here  the 
incident  and  the  reflected  waves  interfere  and  cause  a  periodic 


variation  in  excess-pressure  p  as  well  as  in  p  u2  along  the  axis 
of  the  bean  (see  Eqs.  (2.43)  and  (2.44) s  and  Sec*  II,  7b).  In  a 
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bean  of  infinite  width,  not  in  communication  with  the  undisturbed 
part  of  the  medium  at  the  rear  of  D,  these  tiro  quantities  p  and 


p  u  ,  acting  together  at  the  surface  of  D,  are  of  opposite  sign 
and  cancel  each  other  to  zero  at  small  amplitudes.  But  since 
experimentally  the  beam  is  always  of  finite  width  and  normally 
surrounded  by  the  undisturbed  medium,  the  mean  pressure  p+  inside 
the  beam  is  raised,  as  we  saw  in  Sec.  II,  7b,  up  to  the  static 
pressure  p  ,  Superposed  upon  this  mean  pressure  pQ  is  the 
periodic  variation  p  in  space;  it  amounts  to  p_^  -  pQ  and  is 
different  from  pQ  in  general  at  D.  The  average  negative  value 
of  p  —  p  (which  in  the  infinitely  extended  beam  neutralized  the 
quantity  u2  at  D)  is  now,  in  the  case  of  the  finite  beam,  in¬ 
creased  by  the  amount  E  (l  +  y2)  to  bring  its  value  up  to  pQ; 

p,  -  p  no  longer  neutralizes  the  action  of  the  flux  of  momentum 
t  o 

pu2. 

—  n 

Let  us  write  again  the  expressions  for  p^.  and 
in  a  plane  wave  to  see  what  happens; 


Pt  =  p  +  PQ  =  -  E  (1  +  2  ycos(2kx-+  0)  +  y2)  +  Po  (2.43) 

^>u2  =  +  E  (1  +  2  y  cos(2kx  +  0)  +  (2.44) 

In  the  case  of  a  perfect  absorber  D,  y  =  0  and  there  is  no  varia¬ 
tion  of  p^  and  pu2  with  x.  The  pressure  p^  =  PQ  —  E  is  raised 
to  pQ  in  the  beam  and  all  the  force  on  D  is  caused  by  the  quantity 
p  u2  in  the  way  we  pointed  out  above  under  a.  If  we  consider 
on  the  other  hand  the  case  of  a  perfect  and  rigid  reflector  D, 
there  is  no  movement  of  particles  at  :the  surface  of  Dj  thus  if  the 
origin  x  ~  0  is  taken  at  the  mean  position  of  the  surface  of  D, 
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both  u  and  4  vanish  at  x  =  0,  This  boundary  condition  requires 

that  y-  1  and  9=7 rs  as  is  seen  from  Eq.  (2.24),  from  Eq.  (2.44) 

above  we  find  ^ou^  =  0  and  at  x  =  0.  The  first  terra 

p  u*  is  zero,  because  there  is  no  motion  at  x  =  0$  the  second 

term  p^  =  pQ  because  p  varies  purely  sinusoidally  in  time  at  x  -  0 

around  the  average  value  p  ,  But  there  is  now  the  additional 

effect  of  the  surrounding  medium,  which  raises  the  average  value 

in  space  of  px  to  p  as  illustrated  in  Fig.  14. 
t  o 


2E(1  +  Y2) 


*  %  I  t  / 


e(i  +  y  2)  J-  -  -  -\  -/  -  -  -V  -  -  i~  J 

h,'  i 


Pt  =  P0  +  P 


Pt  =  Po+P  JD 

,  /'  rfi(i  +v2) 

-4 - r — L>  X  ® 


(b)  .. 

Fig.  14 

Time-average  excess  pressure  p(x)  along  the  x-axis  in  a. 
liquid,  when  tire  acoustic  beam  striking  the  device  D 
undergoes  perfect  reflection  at  D  •('£'=r  "!)•  (a)Beam  of 
infinite  width  and  not  in  communication  with  undisturbed 
regions  of  liquid.  The  average  value  in  space  of  pi  is 
lower  than  p0  by  the  amount  E  (l  +  y2)  =  2E.  The  value  of  ' ' 
p  at  x  =  0  is  zero.  (b)Beam  surrounded  by  or  in  communica¬ 
tion  with  undisturbed  regions  of  liquid.  The  average  value 
in  space  of  p^is  raised  by  E  (l  +  =  2E  and_therefore  _ 

equal  to  p  ,  . .  The  value  of  p  at  x  =  0  is  now  E  (1  +  -  2E 
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The  amount  by  which  p^  is  raised,  to  bring  its  average 
value  in  space  to  p^,  is  +  E  (l  +  y2),  as  we  found  already  in  Sec. 
II,  7b.  At  the  front  of  D  appears  now  a  mean,  pressure  (p^)^  — 

P0  +  E  (1  +  y2) ,  or,  with  y  =  1  for  a  perfect  reflector,  (p^_)  = 

pc  +  2  E.  At  the  rear  of  D,  the  pressure  is  pQ.  Consequently 
the  force  exerted  per  unit  area  of  the  receiving  part  of  D  equals 
2  E,  as  expected.  In  the  case  of  a  perfect  and  rigid  reflector 
the  radiation  pressure  measured  experimentally  is  due  therefore 
to  a  mean  excess  pressure  at  the  receiving  plane  rather  than  to 


the  quantity  ^su2,  which  is  aero. 

c.  In  the  general  case,  where  the  incident  energy  is 
partially  absorbed  and  partially  reflected,  the  force  exerted  upon 


D  is  due  to  both  effects;  the  quantity  ^>u2  and  the  mean  hydro¬ 
static  excess  pressure  p  —  p  .  Their  sum  always  amounts  to 

Tr  O 

E  (1  +  \j2),  but  the  amount  contributed  by  each  component  depends 
on  the  reflection  coefficient  y  and  the  phase  angle  9  of  the 
reflected  wave.  From  Eqs.  (2.43)  and  (2.44)  above  we  find  the 
following  values  for  P  =  p  +  E  (l  t^2)  at  E,  where 

x  =  0: 


"x=0  =  -  2  y  E  cos  9  +  pc 


(2. 


or 


=  +  2  y'E  cos  9  +  E  (1  +y2) 

P4o  =  5  (T+y)2-4^Esin2  | 


(2. 


For  the  two  special  cases  first  considered  above,  we  find  here 
again  the  previous  results!  At  a  perfect  absorber  with  y  =  0 


2  •— 

we  get  P^q  =  pQ  and  ^u^_0  =  E,  the  whole  force  being  due  to 
flux  of  momentum  pu^_0.  At  a  perfect  rigid  reflector  with 
y  -  1  and  9=7 r  we  obtain  P^g  =  pQ  +  2  E  and  ^)u^_q  =  Os 
the  whole  force  being  due  to  the  mean  excess  pressure  P 


x=0 


The 


foregoing  equations  give  the  contributions  of  P^q  and 
f  ux=0  ’blle  res^ltant  force  pQ  +  E  (H  +  \^)  per  unit  area  of 
the  receiver  for  all  values  of  \/  and  9. 

Thus s  in  the  general  case*  where  a  reflected  wave  exists 9 
one  has  strictly  to  recognize  these  two  different  effects  acting 
simultaneously  at  a  s unface  struck  by  a  compressions!  wave. 

What  is  measured  experimentally  as  radiation  pressure  is 
the  sum  of  the  flux  of  momentum  and  the  excess  pressure  at  the 
receiving  area,  which  amounts  always  to  E  (1  +  y^)c  The  two  terms 
of  the  sum  depend  on  the  reflection  coefficient  and  the  phase 
angle  9  of  the  reflected  wave 5  but  their  sum  is  independent  of  9 
and  simply  equals  the  sum  of  the  incident  energy  density  E  and 
the  energy  density  E  of  the  reflected  wave.  This  justifies* 
so  far  as  numerical  results  are  concerned;,  the  frequently  adopted 
point  of  view*  which  regards  the  incident  wave  energy  as  perfectly 
absorbed  by  the  receiver*  leading  to  a  contribution  E*  and  then 
partially  or  totally  reenitted  by  the  receiver*  leading  to  an  ad¬ 
ditional  force  ^  ^  E,  This  concept  gives  the  right  numerical  value 
for  the  radiation  pressure^  it  does  not  provide  however  the  physical 
background  for  the  forces  really  acting  at  the  receiver.  (See  also 
p.  49) 

In  the  case  where  the  acoustic  radiation  traverses  a 


plane  reflector  separating  two  different  media  1  and  2*  some  of  the 


energy  being  transmitted  into  medium  2  behind  the  reflector, 
the  resultant  radiation  pressure  acting  upon  the  reflector  at 
normal  incidence  is  obviously  given  by 
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?rad  pt1  ~  pt2  +  P  1ui  ”  f  2U2 


(2.78) 


If  the  media  1  and  2  are  under  the  same  hydrostatic  pressure  pQ, 


then  p.^  —  p^2  =  P-j  -  Pq.  Per  small  amplitudes  and  a  beam  of 
finite  width,  we  can  write  p1  +  p^u^  =  (1  +  ^  2);  assuming 

that  no  reflected  wave  exists  in  medium  2,  we  may  introduce 
the  amplitude  transmission  coefficient  S  by  E0  =  P  u2  =  & 


C2  P  2V 


and  obtain 


P  =  E,  (1  + 
rad  1 


'2  -  &  2) 


(2.79) 


as  already  stated  on  page  53.  If  the  reflector  does  not  absorb 
any  energy,  (1  -  y^'2)  c-jE-j  -  O2o^  or  =  —  (i  -  ^  2)  and 
we  obtain  from  (2.79) 


'  rad  1 


(1  -71)  +V  (1  + 


(2.80) 


This  shows  that  the  direction  -of  the  radiation  pressure  can  also 


2  c* 

reverse  its  sign  if  V  — L— — £  and  c^c  •  *  No  resultant 

5  C1  +  °2  1  2 


p  C-j  —  Cp 

radiation  pressure  would  be  found  if  \/  =  - — rp  -Q—  , 

1  2 


If  the 


two  media  are  substantially  the  same,  c^  =  Cg,  and  from  Eq.  (2.30) 


we  have 


P  ,  =  2  V2  E 
rad  ft 


(2.811) 


*  See,  for  example,  the  experiments  of  Hertz  and  Mende' 
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PART  III 


RADIATION  PRESSURE  IN  GASES 


1 .  The  Lagrangian  Nave  Equation  for  Gases < 


Absorption 


The  Lagrangian  wave  equation,  which  is  linear  in  the  case 
of  liquids  due  to  the  assumption  of  constant  compressibility,  be¬ 
comes  non-linear  for  gases ,  as  the  adiabatic  relation  between  total 

V" 

pressure  p*  and  density  p&  has  the  well-known  form  p.^*  V  4  c  =  constant, 


V  1 
vo  ! 


(3.1) 


where  V"  is  the  ratio  of  the  specific  heats.  (This  notation  is 
5  c 

chosen  in  order  to  avoid  confusion  with  the  coefficient  of  reflec¬ 
tion  bp.  Inserting  Eq.  (3.1)  in  the  Lagrangian  wave  equation 
(2.19) s  considering  T*  =  -  pf  and  Eq.  (2.28)  we  get,  disregarding 
viscous  absorption 


Ik  =  .  i- 

3t2  a?  {1  +  +  Yc 

dx 


(3.2) 


with  cr 


Vc' /  p\ 


This  equation  cannot  be  solved  exactly  by  known  functions 
in  a  closed  form.  So  one  usually  is  content  to  treat  only  the 
case  of  small  amplitudes,  using  the  Taylor  development  of  the  term 


(1  +  +  ^  cad 


and  retaining  first  order  terms.  in  dealing 


with  average  time-values,  however,  one  must  in  general  be  careful 

p 

not  to  omit  second  order  terms.  The  flux;  of  momentum  pu  ,  as 


si 


we  have  learned,  involves  no  second  order  terms  of  Eq.  (3*2) 
at  snail  amplitudes  according  to  our  considerations  leading  to 
Eq.  (2.  37).  On  the  other  hand,  the  correct  calculation  of  the 
"Rayleigh"  pressure  p*  requires  the  solution  of  (3.2)  up  to  second 
order  terns, 

2.  The  Radiation  Pressure  upon  a  Perfect  Absorber 
in  Gases  at  Small  Amplitudes  due  to  a  Bean  of 
Finite  Uidth 

Developing  the  right  side  of  Eq.  (3.2)  in  a  Taylor 
series  and  retaining  only  second  order  terms,  one  finds: 


-°2(1+yc>H  H  (3.3) 

Airy  has  given  the  solution  of  Eq.  (3-3)  for  small  amplitudes 
(l<4<0)  iii  the  form  of  a  series  in  powers  of  k£Q  for  a  pure 
sinusoidal  excitation  assumed  to  be  located  at  the  origin,  where 
x  =  Oj  retaining  second  order  terms,  it  is  known  to  have  the 
following  form  for  a  wave  traveling  in  the  positive  x-direction 
(see  ref.  9,  p.  I30)s 


5  =  K  sin  (cot-kx)  +  - - Y0)..*  (],£  )2  r-j  +  cos  2  (oit-kx)  j  (3.4) 

o  g  O  u.  —I 

For  computing  the  radiation  pressure  of  a  beam  of 
finite  width  upon  a  perfect  absorber,  we  have  from  Part  II,  7a, 


P  ,  =  pu^ 

rad  \ 

Only  the  first-order  term  in  Eq.  (3 .4)  is  needed  for  computing 


(3.5) 


Prad  a_t  sna11  ^Plii^es*  To  this  approximation  ^u^  becomes 
identical  with  ^u-2  and  we  find  easily  from  Eq.  (3.4) 


32 


(3.6) 


At  snail  amplitudes  in  cases  the  mean  kinetic  energy  density 
equals  the  mean  potential  energy  density**  Therefore 
2  =  2  E  ^  =  E  =  mean  total  energy  density.  Thus  we 

obtain  at  small  amplitudes  in  gases  the  same  result  that  we 
fotmd  in  liquids; 


(3.7) 


The  radiation  -pressure  upon  a  perfect  absorber  in  gases  for  an 

acoustic  bean  of  finite  width  equals  the  energy  density  for  snail 

amplitudes.  For  larger  amplitudes  we  have  to  use  Eq.  (3.5). 

The  solution  for  finite  amplitudes  in  gases  is  rather  complicated , 

as  it  requires  a  solution  of  the  non-linear  differential  equa- 

7 

tion  (3.2),  which  can  only  be  given  by  a  series  development. 

Its  consideration  is  beyond  the  scope  of  the  present  paper. 

3 •  The  Fayleinh  'Pressure  in  Gases  for  Progressive 
haves  and  Small  Amplitudes 


For  computing  the  mean  value  p*  which  would  be  observed 
at  a  moving  particle,  we  start  from  Eq,  (3«l)  and  obtain  by 
replacing  (p'VpQ)^C  according  to  Eq.  (2.23) 


(3.8) 


*  See  for  instance  ref.  12,  p.  223 


Using  the  solution  Eq.  (3.4)  for  Z  in  the  abbreviated  form 


K  =  (k?Q)  +  Ji2°  (k%0)2  +  .... 


we  get 


n-M1-  To  ^ * «.)*  (0%  g) 1 

(3.10) 

From  Eq.  (3.10)  it  is  evident  that  we  need  the  solution  for  Z 

up  to  the  second  order  terns,  which  involves  the  tern  in  , 

if  we  wish  to  compute  p*  correctly.  It  is  the  term  VI  *  c>%/d  X 

t  GC  * 

which  was  omitted  by  Rayleigh. 

Inserting  and  which  are  given  immediately  by 
Eq.  (3.4)  s  into  Eq.  (3.10),  and  talcing  the  time-average  value  for 


p*,  we  get 
t 


p~  =  p 
t  o 


1  +  0<£o)  y 


2  f  YcS1J2  _  thb  +  Ye)  ;  j 


4 


3 


/ 


(3.11;) 


whence 


Yc(1i  +Yc) 


p?  =  p  — —  <&0r  +  pc 


(3.  IK) 


Introducing  ^  Pq/c^,  k  =  ci)/c,  and  the  mean  total  energy 

density  E  =  2.  o?  we  obtain  for  the  Rayleigh  pressure* 

in  a  progressive  wave 


4  E 


(3.13) 


This  formula  differs  from  Rayleigh* s  well-known  result  (Hc  +  E/2. 


*  This  is  in  accordance  with  results  obtained  also  by  P.  J, 
llestervelt16  and  J.  S.  Mendousse"1  '  |  both  authors  however  did 
not  proceed  to  the  final  expression  for  ___  ,  as  given  in 

Eq,  (3.13)  here.  t  0 
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The  reason  is  clearly  seen  from  Eq.  ( 3 ® 1 1 ) »  in  which  the  term 
—  y*c(l  "**  y^)/^  as  ^ue  ^ie  -^ac^  ‘t,’ie  second-'order  term  in 

Eq.  (3.4)  is  taken  into  account.  Thus  Rayleigh's  formula  must 
be  multiplied  by  the  factor  l/2,  to  be  correct  in  the  case  of 
progressive  -waves.  The  remark  may  be  added  that  the  formula 
(3.13)  holds  only  at  not  too  large  distances  x  from  the  pure 
sinusoidal  source  of  excitation.  The  second  term  of  the  solution 
(3.4)  contains  the  factor  x.  This  corresponds  to  the  known  fact 
that  a  pure  sinusoidal  wave  in  gases  is  distorted  with  increasing 
distance  from  the  source,  the  wave  energy  being  transferred  more 
and  more  from  the  original  frequency  into  higher  harmonics.  Only 
in  the  immediate  neighborhood  of  x  =  0  is  it  permissible  to  define 
energy  density  by  the  expression  we  used  for  E  and  therefore  to 
attribute  a  reasonable  physical  meaning  to  the  relation  (3.13)  as 
well  as  to  the  result  in  Eq,  (3.7)  for  the  radiation  pressure. 

4*  The  llean  Pressure  p  in  Gases  for  Progressive  Waves 
and  Small  Amplitudes 

Under  the  sane  limitations  just  mentioned,  we  can  give 
the  expression  for  the  mean  pressure  p^  in  Eulerian  coordinates 
in  a  gas  traversed  by  a  plane  infinitely  extended  conpressional 
wave.  Prom  the  generally  valid  Eq.  (2.75)  we  have  in  a  progres¬ 
sive  wave  p  =  p*  —  E.  Introducing  Eq.  (3.13)  and  adding  the 
hydrostatic  pressure  we  find* 


or 


P*  =  P„  + 


t  o  4 


E  -  E 


Po.  -  P 


3  rYc 


to  4 


E 


(3. 


*  See  also  P.  J.  Uestervelt,  ref.  16,  Eq.  (29) 


T2ius,  in  gases  ,  the  terras  p  and  p  do  not  cancel  each  other  in 
the  tensor  component  +  pu^,  as  they  do  in  a  liquid  with 
constant  compressibility.  In  gases  we  have  in  the  absence  of  a 


reflected  wave 


— 5  i  +yc  - 

=  pt  +  f  u2  =  p0  +  ~T~  E 


(3.15) 


5.  The  Radiation  Pressure  at  a  Perfect  Absorber  in 


Gases  for  a  Plane 


xtended  Heave  act  Small 


Amplitudes 

For  an  infinitely  extended  plane  wave  the  radiation 
pressure  is  given  by  =  p  +  pu2,  which  at  small  amplitudes 

equals  p  +  E.  Since  .for  a  progressive  wave  p  =  p*  —  E  according 
to  Eq«  (2,75),  ire  find  in  this  case 


-  _-r  -  H+Vc  - 

prad  ~  PRayleigh  4  E 


(3.16) 


The  Rayleigh  pressure  (in  the  corrected  form  for  pro¬ 


gressive  waves)  becomes  therefore  identical  with  the  true  radiation 

pressure  in  a  beam  of  infinite  width,  which  is  not  in  communication 

with  a  gas  region  unaffected  by  the  wave  notion. 

It  would  be  very  difficult  to  realise  this  case  by 

experiment,  even  approximately^  we  usually  have  to  deal  with  a  beam 

of  finite  cross  section  and  in  communication  with  an  undisturbed 

gas  region.  But  in  the  latter  case  the  mean  pressure  p^  in  Eq. 

(3.14)  is  raised  to  p  and  the  measured  radiation  pressure  becomes 
_  o 

identical  with  pu2  =  Es  as  shown  in  Eq.  (3.7). 

6,  Mean  Density  in  Progressive  Haves  in  Gases  at 


Small  Amplitudes 


A  brief  treatment  of  this  subject  is  given,  especially 


to  point  out  the  difference  between  liquids  and  gases 


The  time  average  value  of  the  change  in  density 
can  be  derived  from  Eq.  (2.23).  At  small  amplitudes  we  have 

£  =  (i  +  -  n  -  +  (~£)2  +  ... 

p0  ~bx  "&X  “b-X 

The  solution  for  £  up  to  second  order  terms  is  given  by  (3.4) • 
Inserting  this  solution  in  (3»17)  and  retaining  all  second  order 
terms  leads  to  the  expression  for  the  density  ^  in  Lagrangian 
coordinates.  The  change  in  density  A  p  =  p  -  pQ  in  Eulerian 
coordinates  is  found  by  substituting  x  —  £  for  x,  according  to 
Eq.  (2,21).  Doing  so,  developing  the  result  in  powers  of  , 
retaining  second  order  terms,  and  finally  averaging  in  time;,  \re 
find 

A 


(3 


-J-  =  (IcE  )2 

po  fo  8  ° 


(3. 


The  mean  density  in  an  infinitely  extended  plane  progressive  wave 

in  gases  is  lowered  by  the  amount  given  in  (3.13) ,  that  is  by  a 

term  as  low  as  the  second  order  in  kc  .  In  liauids  of  constant 

-  -=0  - 

compressibility,  on  the  other  hand,  we  found  that  the  change  in 
average  density  vanished  in  the  second  order  approximation,  as 
it  is  only  of  fourth  order  in  k^Q.  (See  Eq.  (2.68a),  p.  63). 

The  mean  pressure  p+  in  a  progressive  wave  is  lowered 
by  the  wave  notion  by  the  amount  given  in  Eq.  (3.14).  In  a  beam 
of  finite  width  surrounded  by  an  undisturbed  medium  of  hydrostatic 
pressure  pQ  the  bean  undergoes  a  slight  compression  until  the 
pressure  inside  the  bean  equals  pQ,  as  pointed  out  already  in 
the  case  of  a  liquid  medium  (Part  II,  Sec,  7a).  The  increase  in 
density  Ap'  inside  the  bean  due  to  this  compression  is  given  by 
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=  -£=_e°  =  (P°-.tA£)1/V ° .  1  =  i  &tt 

Po  Po  Po  y0  P0 


whence.,  according  to  Eq.  (3.14)  s  we  find 


^  p'  _  3  -yc 

Po  4 


®  -LUC 

fo°2  8 


(k?  )2 

o 


(3.19) 


The  resultant  change  in  density  in  a  bean  of  finite  width  in  a 
progressive  wave  in  a  gas  at  snail  amplitudes  is  therefore  given 
by  neans  of  Eqs.  (3.18)  and  (3.19)  as 


Aptotal  =  Af+Af 


I  r: 


1  -id  1-  a  Irish  \2 

o  4  ’o' 


f 


(3.20) 


7.  Radiation  Pressure  in  Gases  at  a  Perfect  Reflector 
due  to  a  Finite  Beam  at  Email  Amplitudes 
The  boundary  condition  at  a  perfect  and  rigid  plane 
reflector  is  by  definition  (  =  0  or  us  =  u  =  0,  The  corresponding 
first-order  solution  of  the  wave  equation  (3.2)  or  (3.3) ,  which  we 
denote  by  and  respectively s  is 

j"(l)  =  (sin (cot  -  kx)  -  sin(cot  +  kx)  ) 
o 

^(l)  =  sin  lac  cos  (cot  +  tt) 


or 


and 


utt(l)  =  — S.  =  2cx%  sin  lex  sin  cot 
o  t  o 


(3.21) 

(3.22) 


the  well-known  expressions  for  a  standing  wave;  it  fulfills  the 
boundary  condition  u  =  0  at  any  time,  the  reflector  being,  for  ex¬ 
ample;,  at  x  =  0.  From  (3*22)  we  obtain  the  mean  kinetic  energy 
density,  averaged  in  time  and  space, by 


2  Ekin  =  fo  u*^  =  Po^=  fo^ 


2r  2 


(3.23) 
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This  expression  is  correct  up  to  the  second  order  of  approxima¬ 
tion  as  already  stated  on  Page  35 • 

If  we  apply  the  method  mentioned  by  L.  Brillouin  and 
outlined  on  Page  49  of  this  report  for  computing  the  radiation 
pressure  upon  a  reflector,  ue  have  to  add  the  radiation  pressure 


•  4~- — ~  s  regarded  as  reenitted  by  a  reflector  of  ampli- 

PcP^Kc 

tude  reflection  coefficient  to  the  radiation  pressure  - ~ — 

found  in  Eq.  (3.6)  as  due  to  an  incident  and  perfectly  absorbed 
plane  wave.  The  sura,  results  in 


P 

rad 


a  +  f  >  • 


(3.24) 


and  for  the  perfect  reflector  uith  y  =  1  we  obtain 

d  ,  P ow  ?o 

Y=1 


‘  rad 


'2r  2  -  2  5lcin  =  2  E 


(3.25) 


where  E  denotes  the  mean  total  energy  density  of  each  of  the  two 
progressive  uave  components  leading  to  Eq.  (3*21 ),  in  conformity 
with  the  usage  in  this  report,  according  to  Eq.  (2.42). 

The  problem  of  a  perfect  and  rigid  reflector  was  also 
treated  by  F.  Bopp^.  It  nay  contribute  to  one's  confidence  in 
the  result  expressed  by  Eq.  (3.25),  if  we  trace  also  Bopp's 
reasoning  for  this  case,  especially  as  it  can  be  done  very  briefly. 
Averaging  in  tine  the  generally  valid  Eulerian  equation  (2.14) s 
we  have  for  a  wave  notion  periodic  in  time  "3  (p  +  ^u^/bx  =  0, 
and  integrating  with  respect  to  x,  we  find 


c 


(3.26) 
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At  the  perfect  reflector  u  =  0;  thus  the  constant  C  equals  the 

mean  excess  pressure  p  at  the  reflector,  which  is  just  identical 

with  the  radiation  pressure  (P^)  _  in  our  case.  (See  Part 

1 

II 9  Sec.  12b).  On  the  other  hand  we  obtain  by  averaging 
Eq.  (3.22)  in  space 


p  +  jo^  =  G  =  p  +  2  Ej,in  (3.27) 

For  the  beam  of  finite  width  under  action  of  the  static  pressure  pQ 
of  the  surrounding  medium  it  is  reasonable  to  assume,  in  agreement 
with  Bopp,  that  p  =  0  (See  also  Part  II,  Sec.  7b),  in  other  words 
that  the  average  value  in  time  and  space  of  the  total  pressure  p^ 
inside  the  beam  equals  the  static  pressure  pQ  outside  the  beam. 

With  p  =  0  we  get  from  (3.27)  C  =  (Pracj)  =  2  in  agree- 

ment  with  Eq.  (3. 25).  (See  also  Page  92*Eq»  (3.41)) 


8.  The  P-adlation  Pressure  in  Gases  at  a  Perfect 

Reflector  Caused  by  an  Infinitely  Extended  Plane 
Wave  at  Small  Amplitudes 

The  solution  of  this  problem  encounters  some  difficulties, 
as  it  turns  out  that  the  complete  solution  for  a  standing  wave  in 
gases  cannot  be  achieved  without  taking  into  account  the 
dissipation  of  energy,  which  actually  is  of  course  always  present. 
Nevertheless  it  is  possible  to  compute  the  radiation  pressure  for 
this  case  without  proceeding  to  the  complete  solution,  as  will 
now  be  shown.* 

The  radiation  pressure  at  the  perfect  reflector  becomes 
identical  with 


(5rady=  i  =  -  Pc  =  p*  =  p)^ 


(3.28) 


*  See  also  Westervelt's  paper,  ref.  16,  Sec.  ¥1 


90 


since  all  these  pressures  are  the  sane  at  the  inmovable  surface 

of  the  reflector.  The  tine  average  value  of  p*  is  given  by 

t 

Eq»  (3. 1 0)  i 


At  _ 

p£-po 


From  the  first  order  solution  in  (3.21 )  we  find  easily 


(~ — -)=  0  and  — 1)  =  1  +  cos  2  kx  (3. 30) 

3x  c>x  / 

1 

To  find  the  tern  ~£>J?  /<hx,  which  is  related  to  the 

-  2 

second  order  tem  H^2)  of  £  =  ^(^)  +  £'l2^,  we  need  the  solution 
of  the  differential  equation  (3.3) s  which  in  this  case,  introducing 
the  first  order  solution  into  the  last  term  on  the  right  side, 

becomes 


A^OL  =  c2  -£~Ja  +  2  c2  (1  +Y*  )  k^£  2  sin  2  lex  cos2  cot  (3.31) 

^t2  "bx2  «cy  so 

A  solution  compatible  with  the  boundary  condition 
^  =  0  at  x  =  0  at  any  t,  is  found  to  have  the  form^ 


2"  =  +  ^2)  _  j-(l)  +  f(t)  .  (k?Q ) 2  skn  2  lex  (3.32) 

Inserting  (3*32)  into  (3.31)  with  given  by  (3.21),  xre  obtain 
for  the  time-dependent  function  f(t)  the  differential  equation 

-2 

+  4  oi2f  =  (-]  +y  )  coc  (1  +  cos  2  cot)  (3.33) 

dt2  a  c 

This  is  the  differential  equation  of  a  resonant  system  driven  by 
an  impressed  force  of  exactly  the  resonance  frequency  2tn.  As 
there  is  no  dissipation  term  involved  on  the  left  side  of  (3.33), 
the  final  amplitude  becomes  infinite,  as  is  well  known.  This 
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fact  can  easily  be  understood;  Due  to  the  non-linear  relation 
between  pressure  and  density  in  gases,  the  fundamental  angular 
frequency  co,  which  excites  the  standing  wave ,  sets  up  higher 
harmonics  of  angular  frequencies  2co,  3co,  ...  in  the  system.  The 
amplitudes  of  all  the  harmonics  are  limited  actually  by  the 
unavoidable  dissipation  of  energy  in  the  resonating  system,  which 
hov/ever  is  not  taken  into  account  in  the  differential  equation 
(3.31).  A  solution  leading  to  finite  amplitudes  can  only  be 
achieved  if  dissipation  or  absorption  of  energy  is  properly  taken 
into  consideration. 

In  our  problem,  however,  we  are  only  interested  in  p*, 
which  involves  tine-average  values  of  Fortunately  the 

time-average  values  of  all  those  purely  sinusoidal  higher  harmonics 
vanish;  the  only  term  of  f(t)  in  the  differential  equation  (3.33) 
that  is  of  importance  for  us  here  is  the  constant  term  (1  +  ^)  coc 
on  the  right  side  of  (3.33),  which  leads  to  a  time-independent  term 

f  =  ^  =  1  +Yc  (3.34) 

4  co  4  k 

Uith  the  aid  of  (3.32)  we  have  therefore 


?(2)  =  (k?0)2  i2(x)  =  (k^)2 


L sin  2  kx 
4  k 


and  hence 


From  Eqs. 


(3.30) 


'ciJl  2  1  + 

c!x 

,  (3.36)  and 


— cos  2  lex 
2 

(3.29)  we  obtain 


(k^0)2  (l  +  cos  2  lex  -  cos  2  kx) 


or 


£ 


0^0  2 


(3.35) 

(3.36) 


(3.3?) 
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Talcing  account  of  (3.23)  and  (3.28),  ue  have  finally 


'had> 


=  =  (1 +yc)  lcln  =  (1 +^c)  s 


(3.38) 


It  is  interesting  to  note  that  the  "Rayleigh  pressure" 
as  observed  at  a  noving  particle,  becones  independent  of  x 
in  a  standing  uave. 

From  (3.37)  ue  are  able  to  derive  the  time-average  value 
of  the  excess  pressure  p  in  Eulerian  coordinates*  uith  tlie  aid  of 
Eqs.  (2.74)  and  (3.1 8)  ue  find 

_  _  i  +  \/ 

p  =  p*  -  2  E]cin(x)  =  ( . ~'~2' ~ ~  ~  1  +  cos  2  lex) 

or  p=E  (y  -  1  +  2  cos  2  kx)  (3.39) 

kin  ft  c 


The  average  value  in  tine  and  space  of  p  is  therefore 


p  =  (Yf 


1)  E. 


kin 


(3.40) 


The  mean  change  in  pressure  due  to  the  uave  motion  is  positive  in 
gases ?  since  ^  7 1 .  If  the  beam  is  surrounded  by  an  undisturbed 

medium  uith  hydrostatic  pressure  pQ,  it  tends  to  expand  a  little, 
until  the  average  total  pressure  p^equals  the  hydrostatic  pressure 
pQ  outside  the  beam.  In  the  beam  of  finite  uidth  the  pressure  p 
becomes  therefore  according  to  (3. 39) 


p  =  2  E,  .  cos  2  lex 
*  kin 


(3.41) 


At  the  reflector,  xjliere  x  =  0,  ue  have  then  (p)  „  =  2  E_  .  ;  this 

x=0  kin 

leads  us  back,  by  a  different  route,  to  the  result  found  in  (3.25) 
for  the  radiation  pressure  in  a  beam  of  finite  uidth. 
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9.  Note  Concerning  Rayleigh's  Origins!  Formula 
on  the  "Pressure  of  Vibrations" 

Rayleigh’s  original  formula  (energy  density)1-^ 

is  found  in  many  books  and  publications  as  the  formula  for  radia¬ 
tion  pressure  in  gases,  Nevertheless  it  does  not  express  the 
radiation  pressure  as  usually  measured  in  a  beam  of  finite  width, 
but  is  rather  associated  with  the  mean  excess  pressure  p**  =  p*  —  p^ 

observed  at  a  moving  particle,  which  may  properly  be  called  the 
"Rayleigh  pressure".  On  the  other  hand  we  found  in  Sec.  3  and  $ 
above  the  following  results  for  the  radiation  pressure  in  an  in¬ 
finitely  extended  plane  wave  in  gases  at  small  amplitudes: 


(P_  J  =  1-iXl  E  = 


'rad 


0  4 


4 


2r  2 
o 


(3.13) 


PadL  ,  -  0  +*,>  E  '  (1  ■%>  P 


^2K n 2 


(3.38) 


T 


One  must  be  careful  as  to  the  meaning  of  energy  densityj 

throughout  this  report  E  denotes  the  mean  energy  density  of  a 

progressive  wave.  In  the  formula  (3.13)  E  becomes  identical  with 

the  actual  total  energy  density  for  a  progressive  wave.  In  the 

formula  (3.38)  the  actual  total  energy  density  equals  2  E  because 

the  energy  density  of  the  reflected  wave  adds  to  the  energy  density 

of  the  wave  incident  on  the  perfect  reflector.  If  ire  introduce 

the  concept  of  the  total  energy  density  E  ...  the  two  formulas 

- _  total 


may  be  written  in  the  form 


^Prad' 


Y=0 


1  +V„  = 

_  A  C  Tp 


Jtotal 


^rad)y  _ 


1  + 


^2  E. 


total 


(3.13) 


(3.38) 
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This  shows  that  Rayleigh's  formula  in  its  original  Tom  applies 

indeed  to  the  radiation  pressure  on  a  perfect  and  rigid  reflector 

in  a  gas  exerted  by  an  infinitely  extended  plane  wave,  which  does 

not  communicate  with  a  medium  under  a  constant  hydrostatic 

pressure  p0  if  the  total  mean  energy  density  E+n^.a-|  means  the 

xdiole  energy  density  E,  .  +  E  =2E  in  the  standing  wave. 

icin  pot  kin 

For  a  progressive  wave  falling  upon  a  perfect  absorber ,  however , 

the  correct  value  is  half  as  large,  if  by  "total  mean  energy  density”, 

1  ,  one  understands  the  total  energy  density  E]c±n+  E  Ejcin+  Epot 

total 

of  the  progressive  wave# 


asss:ely  of  chief  equations  for  liquids  am?  gases* 

A.  General  Expressions  Valid  at  Large  or  Small  Amplitudes 
in  Li quids  and  Gases  (Disregarding  Viscosity) 

Uhen  a  plane  acoustic  wave  of  finite  or  infinite  width  in 
a  medium  (l)  falls  normally  on  a  plane  slab  of  any  material  or  thick¬ 
ness,  separating  medium  (1 )  from  some  other  medium  (2),  the  resultant 


time-average  pressure  upon  the  slab  is  given  by 

5  =  ^  "  Hz +  "fm?  -  f?? 

Medium  (2)  may  be  bounded  by  still  another  medium  which  causes  a 
reflected  wave. 

If  the  two  media  (1 )  and  (2)  are  both  under  the  same 
hydrostatic  pressure  pQ  when  undisturbed,  Eq.  (2*7 3)  reduces  to 


(2.78) 


*  For  Symbols  see  pages  ii  and  iii 
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P  ,  =  p  ~  p  +  £  U  *  «.  p  u 
rad  1  2  fl  1  12  2 


(2. 78a) 


which  now  represents  radiation  pressure  alone. 

B .  Ho  Reflected  Have  Present  in  Iledium  (2) 

1.  If  no  reflected  wave  is  present  in  medium  (2),  ’out 
the  slab  still  transmits  and  -absorbs,  tre  have  for  a  beam  of  finite 
width  at  snail  amplitudes  and  surrounded  by  undisturbed  regions  (l ) 
and  (2)9 

P  =  (1  +  V2  -  6 2)  £  (2.79) 

rad  0  1 

2.  If  no  energy  is  absorbed  within  the  slab,  the  fore¬ 
going  formula  for  a  bean  of  finite  width  at  snail  amplitudes  can 
be  reduced  to 


'  rad 


(1  -  SI)  +  V2  (1  +  21) 

°2  3  c2 


E 


1 


(2.80) 


G.  Opaque  Slab-  Ho  Have  I  lotion  Present  in  Medium  (2) 

1.  If  the  slab  absorbs  all  energy  that  is  not  reflected, 
so  that  no  wave  notion  is  transmitted  into  nedium  (2),  the  radiation 
pressure  exerted  by  a  bean  of  finite  width  at  small  amplitudes  in 
liquids  and  gases  is  given  by 

P  ,  »  (1  +  V5)  1,  (Eqs.  (2.50),  (2.54/9  (2.55), 
rad  5  1  (3.7),  (3.24)9  (3.25)) 

2.  Under  the  sane  circumstances,  but  if  the  acoustic 
plane  wave  is  regarded  as  inf ini t elv  extended  and  not  communicating 
with  undisturbed  regions,  we  have  for 

a.  Liquids  (with  constant  compressibility) 

i?rad  =  0  (Pages  33,  36,  45) 

b.  Gases  (for  adiabatic  processes) 


^rad  =  ^  a  wave  (3*13) 
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(3.33) 


^rad  ~  ^  +y<3^  in  a  standing  wave 


The  symbol  E  stands  for  the  tine  average  of  the  energy  density 
of  a  purely  progressive  acoustic  wave  of  angular  frequency  co  and 


maximal  amplitude  £  in  a  medium  of  undisturbed  density 


Pc 


as 


given  by 


i,  =  .  £ 

1  2  c 


(2.42) 


where  J  denotes  the  intensity  of  the  progressive  wave. 

D.  I  lean  Excess  Pressure  Due  to  a  Plane  Infinitely 

Extended  Acoustic  have  at  Small  Amplitudes 

The  hydrostatic  pressure  p  in  the  undisturbed  medium 

is  changed  by  the  acoustic  wave  motion  into  p  +  p,  where  p  is 

o 

called  the  excess  pressure* 

1 •  Liquids  (with  constant  compressibility) 

p  =  -  i  (1  +  2  Y cos (2  lex  +  9)  +  y2)  (2.43) 

P  =  -  E  (i  +  y2) 

2.  Gases  (for  adiabatic  processes) 


3-y0_ 

— 1  ^  '71 


P  =  ~ 


S  in  the  progressive  wave  ('y'  =  0) 


P  -  L  ( Yl  ~1  +-2  cos  2  las)-  in ’the  standing 
®  wave  (y=  1;) 

p  =  E  ( -  1 )  in  the  standing  wave 


(3.14) 

(3.39) 

(3.40) 
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